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I. INTRODUCTION AND SUMMARY 


Free electron lasers have many attributes that makes them attractive 
for a variety of applications. The attributes include wavelength agility, 
high efficiency and scalability to high powers. The short-wavelength free 
electron laser device may be configurated in two major ways which are char¬ 
acterized by the type of electron accelerator used: radio frequency (RF) 
or induction linear accelerator. Both of these classes of devices have 
optical problems that are unique and, when approached using conventional 
techniques, are technologically difficult. The objective of this program 
was to examine nonlinear optical techniques for addressing these problems. 

This report begins with a summary of the operational characteristics 
of induction and radio frequency driven free electron lasers and conven¬ 
tional approaches to the optical issues. A discussion in more details is 
given again in Section II. The following section describes various non¬ 
linear optical process and delineates areas in which these can be applied. 
As a result of the examination of technical issues and applicability of 
nonlinear optics to these issues, three major areas were pursued. These 
techniques include the multiple-pulse stimulated Brillouin scattering, 
short-pulse four wave mixing and the four wave stimulated Brillouin scat- 
teri 1 ''- The concepts developed using these techniques are phase conjugate 
resonator, phase conjugate outcoupling and phase conjugate beam conversion. 

Theoretical models and laboratory experiments were performed for 
initial analyses and verifications of the derived concepts and the asso¬ 
ciated nonlinear optics technology. These are discussed in section III and 
IV of this report. As a result of this effort, it was found that the non¬ 
linear optics technology is promising for free electron laser system appli¬ 
cations. Some major optical issues associated with the conventional optics 
can be resolved using the identified nonlinear optics techniques, espe¬ 
cially in the areas of resonator mode stabilization and atmospheric beam 
transmission. 
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1.1 FREE ELECTRON LASER SYSTEMS 


Free Electron Lasers (FELs) have been studied for over a decade since 
the first observation of radiation at Stanford University in 1972. Experi¬ 
ments have progressed from the observation of oscillation with a constant 
wigglerl»2 to the verification of increased gain and efficiency using a 
tapered wiggler.3/4 A key feature of FEL resonators is their extremely 
high saturation flux in the interaction region. This leads to a number of 
presently unresolved issues when high power scaling is envisioned.5-11 in 
particular, scaling analyses show that advances in the state of the art of 
resonator design will be required to achieve high output powers.12-14 This 
program had the goal of performing initial measurements and verifications 
of approaches to deal with the issues associated with the high flux as well 
as the atmospheric propagation of high power FEL beams. Th order to dis¬ 
cuss these issues it is necessary to have an understanding of the compon¬ 
ents of an FEL directed energy system. The system is divided into three 
major subsystems: laser device, transfer optics, and beam director. 

1.1.1 Laser 

RF Approach 

The RF linac FEL is shown schematically in Figure 1. An electron beam 
is produced in an injector by thermionic emission and initially accelerated 
to an energy of around 10 MeV. RF cavities then accelerate the beam to 
approximately of 100 MeV which is required for visible FEL operation. The 
accelerator produces a series of short micropulses 10 ns long separated by 
a small integral number of RF periods. Peak currents are in the several 
hundred ampere range, and at this low current and associated low gain it is 
necessary to operate the system as an oscillator. The beam passes through 
a several meter long wiggler where sinusoidal magnetic fields convert some 
fraction (lO 0 .) of the electron energy to light. The spent electron beam is 
sent to a beam dump or decelerated through a separate set of cavities to 
recover some of the energy. 

In oscillator operation the FEL produces an intense beam of small 
divergence, so that high fluxes impinge on the mirrors. At high average 
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powers It is necessary to utilize some advanced technique to handle such 
fluxes, e.g., very long cavities with exacting alignment requirements, 
intracavity aerolenses to magnify the beam spreading, or grazing incidence 
mirrors to spread the beam over a larger area. Generally, with these 
optical approaches it is difficult to achieve good beam quality and state- 
of-the-art alignment tolerances. In addition methods for outcoupling are 
difficult to devise and generally require the complementation of high power 
gratings. 

Typical requirements for such systems are shown in Table 1. 

Table 1. Typical RF FEL Optical System Requirements 



In summarizing the issues associated with such systems, the following 
are seen to be the key issues: 

• Optical damage threshold 

• Alignment requirements 

• Figure requirements 

• Outcoupling methods 

• Use of intracavity with aerolenses 

• Bandwidth control 

• Production of harmonics and sidebands in the FEL 
Induction [ inac 

A schematic of an induction linac system is illustrated in Figure 2. 

An electron beam of several thousand amps and 50 ns duration is produced in 
the injector by thermionic or field emission cathodes. This beam is 
01-127-86 4 
































accelerated to 300 MeV by a series of single turn transformers. High 
average power is achieved by repetitively pulsing the system hundreds or 
thousands of times a second. Because the length of the electron pulses are 
less than the oscillator length, oscillator operation is impossible but the 
high currents involved allow the system to have very high gain which per¬ 
mits amplifier operation. A master oscillator pulse, produced by a YAG 
system, for example, at 1 joule is amplified to several thousand joules 
while passing through a 100 m long wiggler. The optical pulse is essenti¬ 
ally trapped by the electron beam much as in an optical fiber so that 
divergence of the optical pulse does not occur during the passage through 
such a long wiggler. This has the unfortunate side effect of introducing 
jitter to the output beam. The peak intensities involved are enormous, 
since the power is delivered at the end of the wiggler over a fraction of a 
square centimeter. The damage mechanism for the optics is the pulse energy 
rather than average power In this case, and distances of several kilometers 
are required to allow the beam to spread sufficiently so the surface break¬ 
down of mirrors does not occur. Aerolenses could possible reduce this to a 
more manageable distance. 

Typical requirements for such systems are shown in Table 2. 

Table 2. Typical Induction FEL Optical System Requirements 


Optical Systems 
Input Power 
Output Power 
Repetition rate 
Wiggler length 
Damage threshold 
Divergence length 


Requirements 
1-10 Joules 
1-10 Kilojoules 
1-10 Kilohertz 
120 meters 
10 J/cm^ 


5 Kilometers 
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In summarizing the issues associated with such systems the following 
are seen to be the key issues: 

• Optical damage threshold 

• Alignment requirements 

• Production of sidebands and harmonics 

• Use of magnifying aerolenses 
1.1.2 Transfer Optics 

The purpose of the transfer optics is to condition the optical beam 
for introduction into the beam director. For both FEL systems this 
involves establishing the proper beam diameter and phase curvature charac¬ 
teristics, eliminating any unwanted frequency components in the light, 
correcting any locally introduced phase aberrations and jitter, and compen¬ 
sating for any atmospheric or beam director introduced aberrations which 
are expected to occur in propagation through the beam director and atmos¬ 
phere. The conventional approach for correcting these propagation errors 
is to sample the errors with a beam from a beacon. An aperture sharing 
element is required to separate the outgoing wave from the Incoming beacon 
beam, a task made more difficult by the requirement that both waves be at 
approximately the same wavelengths. The analysis of the phasefront errors 
is performed using a phasefront sensor, and a deformable mirror is used to 
introduce error compensating phase variations into the beam. To perform 
this analysis requires a large number of calculation channels of very high 
speed as well as very complex cooled deformable mirror surfaces. For suf¬ 
ficient phase fidelity, thousands of channels or more will be required. 
Although some work in this area has been carried out successfully at low 
powers, the technical effort remaining to make a high quality, high power 
system is still very large. 

1.1.3 Uplink 

The uplink system is essentially a very large astronomical telescope 
operating in reverse. Significant additional requirements are placed on it 
by the high flux loads, surface figure requirements, and fast slewing 
required. In the RF case director apertures of a couple of meters are 
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sufficient to provide resolution of the target and to reduce the thermal 
flux in the atmosphere to manageable levels. The induction linac's high 
peak power pulse has an additional constraint imposed by atmospheric Raman 
conversion which limits the peak flux to less than 1 MW/cm 2 . This results 
in a requirement for the beam director to be on the order of 10 m in output 
aperture, which may be beyond manufacturing capability for a single unit 
mirror. If multiple mirror segments are utilized, the alignment and sur¬ 
face tolerance control are issues which require significant technology 
development for resolution. 

1.2 NONLINEAR OPTICS TECHNOLOGY 

The nonlinear optics techniques that are candidates for FEL applications 
are summarized in Table 3. The techniques include stimulated Brillouin scat¬ 
tering (SBS), stimulated Raman scattering (SRS), stimulated Rayleigh-wing 
scattering (SRWS), four wave mixing (FWM), and, for completeness, pseudo 
phase conjugation (PPC). Their operational characteristics are also listed 
for comparison. The SRWS technique is attractive for its short response time 
and small frequency shift that are the desirable property for a RF linac FEL 
resonatorl 5 .I 6 . However, very little of SRWS study has been reported to 
date. The reflectivity observed was low and the conjugation fidelity was 
poor. Pseudo phase conjugation using arrays of retroreflective corner 
cubes^ 7-10 provides near-unity reflectivity, white-light and instant 
response, no frequency shift and no minimum power requirement. However, 
there are two drawbacks for PPC. Due to the finite size of the small corner 
cubes, it can not correct the phase variations that has a characteristic 
periodicity smaller than the individual cell dimension. The conjugation 
fidelity can thus be poor in many cases. Second, the field cannot preserve 
its initial polarization state upon the reflection 20 . The changes of the 
polarization could reduce the radiation intensity in the interaction region 
significantly. In addition, the mirror damage problem is more serious for a 
PPC mirror because of the sharp surface edges. As a consequence, both SRWS 
and PPC were not included in our research major efforts. 

For SBS, SRS and FWM, the potential applications for RT Linac and 
induction linac FELs were identified and summarized in Table 4. In gen¬ 
eral, SBS is a very efficient process if the input pump power is well 
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SRS 

<1 

$1 

0.1 
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10-12 

1 

>200 
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High 

High 
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-- 
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above the threshold. SBS has the advantage of a small frequency shift that 
could be buried under the FEL gain spectrum. While its relatively slow 
response presents initial problems for much shorter RF linac FEL pulses 
(10-100 ps), the TBS reflectivity can be highly enhanced using the concept 
of multiple-pulse SBS 2 1. In this case, the phonon field decays slowly 
enough to serve as a starting acoustic field for the next pulse. A theore¬ 
tical study of multiple-pulse SBS was conducted during this program. 

Another major phase conjugation technique studied, FWM22 ( is not a 
stimulated emission process and, thus, has no power threshold. It gener¬ 
ates the conjugate reflected wave by parametric mixing of the three input 
waves where constructive superposition is in the forward and backward dir¬ 
ections. High reflectivities up to 100 times have been demonstrated in 
liquid media. Its fast response and degenerate reflectivity would be suit¬ 
able for RF linac FEL applications. 

For induction linac FEL applications, a new nonlinear optics process 
was identified and called four-wave stimulated Brillouin scattering 
(FWSBS). This process is essentially the SBS effect arranged in a four 
wave mixing configuration23. it has the combined advantage of both SBS and 
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FWM effects as shown in Table 5. This effect was studied in detail in this 
program. The theory and modeling will be reported in Chapter III and the 
experimental results will be presented in Chapter IV. 

SRS is a optically fast phenomenon with a large vibrational frequency 
shift. It has been considered for beam combination, wavelength shift and 
beam clean-up. Research in this area has been carried out in the Raman 
Beam Combination Program and will not be reported here. 

In general, bandwidth and system efficiency are the two major issues 
of these three techniques for FEL application. The FEL bandwidth is an 
unsettled issue even from the FEL physics point of view. Some theories 
indicated that the FEL is an intrinsically wide band laser source. Other 
theories have argued that the output bandwidth can be extremely narrow if a 
proper filter or a narrow band injection laser is used. In this program, 
however, the issue was left open and the bandwidth problem was addressed 
whenever necessary. 
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Table 5. Comparison of Characteristics Among the Three Nonlinear 
Optics Effects: FWM, SBS, and FWSBS 






















II. CONCEPTUAL DESIGN 


2.1 TASK OBJECTIVE 

The objective of the conceptual design task in this FEL resonator 
development program was to study unconventional adaptive optics concepts in 
the operational regime of FEL characteristics and provide basic resonator 
designs, for use in high power FEL systems. To accomplish this objective, 
the conceptual design effort was divided into three major tasks: 

a. Identify the technical issues of conventional resonator designs 
that are expected to affect critically the performance of a high 
power free electron laser system. 

b. Evaluate the potential of various nonlinear optics techniques for 
free electron laser system applications. 

c. Develop basic resonator designs and identify issues associated 
with the employed nonlinear optics techniques that will be subject 
to the following theoretical and experimental studies. 

2.2 APPROACH 

In order to determine the bases of consideration required for the 
conceptual designs of a high power FEL resonator, a three-step approach 
was followed: First, the performance of high power FEL systems using con¬ 
ventional resonator designs were evaluated and the potential problems 
pertinent to both the RF linac and induction linac FEL systems were identi¬ 
fied. Second, the nonlinear optics techniques available for possible 
applications were reviewed and their characteristics associated with the 
FEL operational parameters were discussed and included in the Introduction 
section. Third, a system analysis of the FEL resonator was performed to 
determine the limit of power losses for a feasible resonator design. These 
steps led to the conceptual designs described in Section 3. 

2.2.1 Conventional FEL Resonator Issues 

The resonator configuration used in all demonstrated FEL systems is a 
linear stable cavity^ - ** (Figure 3). It consists of two spherical mirrors 
with their centers of curvature overlapped along the optical axis near the 
wiggler center. At high powers, the resonator requires a very long cavity 
for the intracavity mode to have a sufficiently large mode size at the 
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Free-Electron Laser Resonators 



Figure 3. Conventional resonator design for high-power FEL systems 
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mirror. This long cavity has problems in its engineering aberration control 
and misalignment sensitivity. A ring resonator configuration with two pairs 
of hyperboloid-paraboloid mirrors (Figure 3) has been suggested^ -7 to shor¬ 
ten the cavity materially. The grazing-incidence hyperboloid mirror 
increases the beam divergence and the paraboloid mirror recollimates the 
beam. The two paraboloid mirrors can be replaced by two spherical mirrors 
to create an inverted ring configuration (Figure 3). The last configuration 
avoids the difficulty of fabricating off-axis paraboloid mirrors. However, 
the aberration effect and the misalignment sensitivity remain the same even 
though the cavity is much shorter^-S. Mirror figure requirements are 
typically one-hundredth of a wavelength and the alignment sensitivity angle 
is near 1 /*rad or less. These requirements push the mirror and auto align¬ 
ment technology to their present limits. 

The best conventional resonator design for a high power FEL system 
might be an inverted ring resonator configuration (Figure 3). This con¬ 
figuration has the complexity of the grazing-incidence mirrors, but reduces 
the engineering difficulty of dealing with an extremely long cavity. How- 
ever, most of the optical issues associated with linear resonator designs 
can still be found in this configuration. Those issues that might affect 
the performance of a high power RF linac FEL or induction linac FEL are 
summarized in Table 6 and discussed in the following. 

The issues can be divided into three groups. The first group is asso¬ 
ciated with the optical resonator. The second group is related to beam 
transfer and propagation. The third relates special FEL systems that 
couple a number of FEL devices. 

A. Resonator Optics 

Mirror damage is the first major problem that came with the 
development of a high power FEL system. Although this is a general issue 
for all high power lasers, the utilization of grazing-incidence mirrors 
near the high flux region in a high power FEL system presents a serious 
problems for mirror cooling and optical stabi1ization technologies. This 
problem can be less serious using nonlinear optical processes. The medium 
used to create the process can be liquid, gas or even plasma. Using these 
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techniques, it is possible have a much higher thermal damage threshold than 
a conventional reflector as well as the capability for recovery after 
damage. 

Higher harmonics of the optical beam are generated in the FEL and grow 
with the fundamental wave. This phenomenon has been observed in previous 
FEL experiments? and found to have serious degradation effect on the mir¬ 
rors?. The suppression of higher harmonics is thus important to ensure the 
survival of the optical resonator in the visible regime. In a nonlinear 
optics process, the interaction efficiency is strongly dependent on the 
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intensity. Since the high harmonics have lower gains than the fundamental, 
SBS phase conjugation is a naturally selective process that will suppress 
the growth of shorter wavelength components. 


•5 

\ 


Conventional outcoupling schemes including grating, scraper and 
transmissive outcoupling (Figure 4) have serious problem. At high power 
levels, transmissive outcoupling scheme can not be used due to thermal 
absorption problem and resulting distortions. Grating technology is cur¬ 
rently not able to sustain high intracavity fluxes on its corrugated 
surfaces. Since high power FEL device generally has a very low single-pass 
gain (< 10 percent), the scraper outcoupling scheme becomes impractical for 
either stable or unstable resonator geometries. The output would be a thin 
annular beam plagued with significant phase fluctuations. An outcoupling 
scheme using nonlinear optics is an attractive solution for this problem 
because it can produce a unobscured output beam and is subject to less 
thermal damage than gratings and scrapers. 

With respect to the end mirrors of the resonator, the beam waist inside 
the wiggler is in the far field region. A small amount of phase nonuni¬ 
formity generated on the mirrors would produce a significant intensity 
fluctuation in the interaction region and thus reduce the gain. In general, 
mirror aberrations must be controlled to within one-hundredth of a wave¬ 
length, which is near the limit of current mirror fabrication technology. 
Phase conjugation using nonlinear optics techniques has been well proven for 
aberration control in a resonator. The resulting benefit could be two-fold: 
reduction in the aberration requirements on FEL mirror surfaces and still 
ensuring good beam quality. 

The problem of misalignment sensitivity of a high power FEL resonator 
has long been recognized 1 1. It is due primarily to the long cavity length 
and narrow interaction region, and is not improved by using of diverging 
mirrors to shorten the cavity 8 . The tilt sensitivity angle of the end 
mirrors is approximately one microradian. The resonator thus requires a 
very sophisticated auto alignment system with an adjustment accuracy in the 
sub-microradian range. Mirror tilt is a special case of aberrations and 
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can be corrected by the optical phase conjugation technique. This dynamic 
correction of the mirror tilt or jitter can reduce the alignment require¬ 
ment significantly. 

B. Up-Link Propagation 

For a ground-based high power laser system, the transport of the output 
laser beam through an optical train and turbulent atmosphere to relay optic 
in space presents a major challenge 1 ^. One of the solutions is to place the 
oscillator in spaced, a beacon signal coming down from the oscillator 
would sample atmospheric turbulence. The signal would then be phase conju¬ 
gated and amplified to high power. The high power beam would be transmitted 
back into space retracing the optical path thereby compensating the aberra¬ 
tions in the optical train as well as in the atmosphere. This simple scheme 
will not be effective for a FEL system because the phase information will be 
completely lost in the high-gain amplification process. A different techni¬ 
que for imprinting the conjugated phase information onto the high power beam 
is required for an up-link system. 

C: Phased Arrays 

As a ground-based laser system, a single induction linac FEL is able 
to produce a very high power beam. However, in the near future, it might 
be necessary to couple several RF linac FEL devices together to reach the 
required power level 1 ^. Therefore, phased FEL arrays are an active 
research topic which includes different schemes such as injection locking, 
coupled cavities and beam combination . Since the phased-array concept has 
been actively pursued in other HEL programs, it is not included in the 
research scope of this project. 

Most of the issues in the resonator and up-link optics are unique to 
the FEL system. Novel concepts have been explored for developing superior 
solutions. Recent developments in nonlinear optics technology have 
demonstrated control of the mirror aberrations and cavity alignment using 
phase conjugation phenomenon 10 . The feasibility of this approach for FEL 
applications has been addressed in this program. 

In a phase conjugate FEL resonator, the two end mirrors are replaced 
by two phase conjugate mirrors (Figure 5), which are created by nonlinear 
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Phase-Conjugated FEL Resonator 



Figure 5. Phase-conjugate REL resonator configuration with 
conjugators at both ends 























optical interactions in a nonlinear medium. The nonlinear effect that is 
most suitable for this application is stimulated Brillouin scattering 
(SBS). Because the pulse length in a RF linac FEL is much shorter than the 
SBS response time (~ nsec), a multiple-pulse SBS approach was investigated 
in performing the feasibility study. 

For a ground-based laser system, the up-link propagation problem 
through a turbulent atmosphere is another major issues. The concepts of 
phase conjugate outcoupling for an RF linac FEL and phase conjugate beam 
conversion for induction linac FEL were developed in this program. The 
nonlinear optics techniques utilized for the two concepts are short-pulse 
four wave mixing and four-wave stimulated Brillouin scattering, respec¬ 
tively. These two effects have just begun to receive extensive attention 
and prove to be promising for the high power beam up-link from ground-based 
FEL systems. 

2.2.2 System Analysis of FEL Resonator 

Consider a FEL resonator design utilizing ideal phase conjugate 
mirrors. In this case a FWM configuration must be ruled out as a candidate 
because it requires two higher-power pump beams and more power is lost 
overall than what is reflected. An attractive candidate might be a SBS 
mirror. However, there are some SBS losses that would impair the resonator 
characteristics. In order to study this issue with regard to its impact on 
the device performance, system analyses were conducted for RF linac FEL and 
induction linac FEL devices. 

From the theory of tapered wiggler FEL^ 4 , the variation of the 
synchronous particle energy can be calculated from the equation 

H? = A ~~ 4 c - sin * (2.3-1) 

where 7 mc^ is the synchronous particle energy, u s is the optical frequency, ♦ 
is the synchronous particle phase, A w and A s are the dimensionless amplitudes 
of the magnetic field and the optical.field. Assume the change of 7 is 

















small. A w and are kept constant. The change rate of 7 is found to be 
proportional to the square root of the optical power, P, 


„ „ pl /2 
dz K 

Since the FEL interaction gain and efficiency are defined as 

r - A£ 

G " r 

- . . 42 

V p ^ 

eb ? * 

the Equation (2.3-2) can be integrated to obtain 


(2.3-2) 


(2.3-3) 


(2.3-4) 


(G + 1 ) 1/2 - 1 = 


(2.3-5) 


where k is a proportional constant. Equation (2.3-5) was derived without 
any restriction on the magnitude of the gain. For an induction linac FEL 
amplifier, the gain is much larger than 1 . Equation (2.3.5) gives the gain 
dependence on power, 


G = r 


(2.3-6) 


and 7 is a constant. For a RF linac FEL amplifier with a gain much lower 
than 1 , the gain-efficiency relation is 


. 4k 
g 7 = p— 
K eb 


(2.3-7) 


Similarly, for a RF linac FEL oscillator where the power is not at the 
resonance, the relation can also be found to be 


G >7 3 ' k ' 4 P eb 


(2.3-8) 


where k' is another constant 

The analysis for FEL amplifiers is straightforward and is not of 
interest here. However, the situation is more involved for a FEL 
oscillator because of the dynamical response of the intracavity power to 
changes in the FEL gain. To complete the equations required for the system 
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analysis, the gain and the system efficiency, 7 S , are expressed in terms of 
the outcoupling fraction, C, and other losses, L, 

G = , { ^ E (2.3-9) 


n - 7T- (2.3-10) 

K eb 

Now, with the help of gain-efficiency relations, the system efficiency 
can be found to be 

V s = k P 1/4 cTT (1 “ C • L)1M (2.3-11) 

The result is shown in Figure 6. The curve for L = 0 is an ideal condition. 
The value of rj s is normalized to have the maximum of 1 at C = 0 on this 
curve. For L > 0, rj s is zero when C = 0orC + L = l. The curves show that 
there is an optimum outcoupling fraction for any given losses. However, no 
significant variation was found between C = 0.2 and 0.6. For the loss of 5 
and 10 percent, it was found the efficiency is reduced by 20 and 30 percent, 
respectively. 

2.3 SYSTEM CONCEPTS 

Based on the problems identified for conventional FEL resonators and 
our understanding of the unconventional adaptive optics techniques, three 
novel resonator concepts have been developed using the nonlinear optics 
technology. The first one is the phase-conjugate resonator concept that 
would correct the aberration effects and misalignment while also allevi¬ 
ating the mirror damage problem. The high reflectivity resonator mirrors 
are created by the multiple-pulse SBS process. The second concept is the 
phase-conjugate outcoupling of a RF linac FEL resonator using the four wave 
mixing process. It simplifies the problem of outcoupling optics and pro¬ 
duces a conjugate beam. The third concept is the phase-conjugate beam 
conversion of the induction linac FEL device using the FWSBS technique. It 
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converts the high power output beam into the conjugate of a probe beam with 
high reflectivity and high efficiency. 

2.3.1 Conjugate PEL Resonator 

The simplest conjugate PEL resonator design replaces both mirrors of 
the linear resonator configuration with two gas medium SBS cells. Since 
gases can sustain much higher fluxes than mirrors, the SBS cells can be 
placed near the wiggler and thus reduce the overall cavity length sub¬ 
stantially. The reflectivity from the SBS cells may be sufficiently high 
to permit oscillation. However, the phase conjugation property that is 
important for aberration control and misalignment corrections, might be 
absent completely in this configuration. In order to a phase conjugate 
wave, the pump beam has to be focused so that the phase information on the 
beam transformed into an intensity distribution in the focal volume- 
interaction region. The reflected field following this intensity variation 
would have the highest SBS gain, thereby generating the conjugate of the 
input beam. 

A concern associated with the use of gaseous medium SBS cells is 
optical breakdown. The optical breakdown in gases consists of two steps, 
initial ionization and avalanche ionization. At an intensity of 100 
Gw/cm 2 , multiphoton ionization process can readily take place in gas mole¬ 
cules. For low pressure gases and short laser pulses, the secondary 
effect of ionization by electron collisions with atoms that causes the 
avalanche, might be negligible. However, a complete breakdown may occur 
for intensities of 1 TW/cm 2 . For intermediate power RF linac FEL devices 
anticipated to be developed in the near future, SBS cells using inertial 
gas media, such as Xenon, may be sufficient for the application. For 
higher power devices, plasma media should be used. However, little is 
understood of the SBS effect in plasma as a phase conjugation process. 

More studies, particularly experimental work should be conducted before a 
conclusion can be drawn. 

Since the phase conjugation process in SBS requires a focused beam, 
concave grazing-incidence mirror can be used to focus the beam into the SBS 
cell (Figure 7). The mirrors are configured as a ellipsoid with foci at 
the wiggler and at the phase conjugate cell. A spherical mirror can be 
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Grazing Incidence Ellipsoidal Mirror for Focusing 
the Radiation into the SBS Cell 
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placed at some distance behind the cell . This spherical mirror functions 
as a conventional resonation mirror at oscillator start-up and at low 
powers when there is no SBS action. As the intracavity power increases 
above the SBS threshold, the optical wave is increasingly reflected at the 
SBS cell. The spherical mirror then intercepts escaping radiation and 
reflects the wave into the SBS cell as a back-injection signal to enhance 
the reflectivity. 

A complete phase-conjugate FEL resonator can be formed by placing the 
set of elliposoidal-SBS-spherical mirrors on the both sides of the wiggler 
(Figure 8). A grating outcoupler is included just to show a complete con¬ 
figuration. In general, use of gratings is less favorable inside high- 
power resonators. 

The pulse length in RF linac FEL devices ranges from 3 to 100 psec, 
depending on the accelerator type. However, the characteristic response 
time of the SBS effect is in the nanosecond regime. The laser field does 
not stay in the SBS medium long enough to permit the Stokes wave to build 
up from noise. However, if the pulse-to-pulse interval is comparable to or 
shorter than the phonon relaxation time, the acoustic wave will persist in 
the medium until the next pulse arrives. A significant SBS reflectivity is 
thus possible even with extremely short pump pulses^l. This multiple-pulse 
SBS phenomenon is essential for a successful conjugate FEL resonator and is 
a subject of theoretical studies in this program. 

2.3.2 Phase Conjugate Outcouplinq 

The identification of a satisfactory FEL outcoupling scheme has been a 
major problem in designing a high power RF linac FEL resonator. A grating 
can produce an unobscurred beam. However, there are two output directions 
and, gratings exhibit much lower damage thresholds. Scraper outcoupling is 
practical only if the outcoupling fraction is significant. For most RF 
linac FEL systems, the gain is approximately 10 percent. The result is an 
undesirably thin annular beam. To solve this problem, a phase conjugate 
outcoupling scheme using nonlinear four wave mixing process has been 
identified. 
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The FWM cell, containing nonlinear gas media, is placed inside the 
resonator (Figure 9). Its position is adjusted such that the counter- 
propagating short pulses overlap inside the cell. There is no stringent 
alignment problem of the pump pulses because only the on-axis pulses can 
survive the oscillation. A low power probe signal with aberrated phases is 
incident into the interaction region. The nonlinear mixing process gener¬ 
ates a conjugated field propagating in the direction opposite to the 
incoming signal. This outcoupling scheme has the following advantages: 

a. The output is a well-filled beam suitable for propagation 

b. The reflectivity can be much larger than 1, such that a weak 
signal can induce a high power return26 

c. Typical maximum energy transfer efficiency of FWM is 5 to 10 
percent which is roughly equal to the FEL single-pass gain 

d. The reflected wave is a conjugate of the signal and can correct 
the phase aberration on its return path 

e. Typical response time of FWM are several picoseconds, which are 
much shorter than the pulse length 

f. The outcoupling is negligible at low powers which can shorten the 
time for the intracavity field to reach saturation 

g. The damage threshold of gas media is much higher than conventional 
mirror surfaces. 

However, there is an equal amount of energy scattered in the forward 
direction accompanied with the conjugate reflection. This reduces the out¬ 
coupling efficiency by 50 percent. Several designs have been considered to 
alleviate the problem. One possible solution is to inject the forward wave 
back into the resonator for further oscillation. The second concept is to 
reflect the forward wave by another phase conjugate mirror, such as an SBS 
cell. Therefore, there are two conjugate return pulses for each interaction 
The third method takes advantages of the short pulse length and finite relax 
ation time of the index grating formed in the medium (Figure 10). By adjus¬ 
ting the arrival time of pump pulses with respect to the signal, it is 
possible to create different gratings to reflect the pump beams selectively 
into the forward or backward direction respectively. It is called 
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Figure 9 


Phase-conjugate outcoupling from a FEL resonator, using the 
four wave mixing effect 
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"Discriminative scrattering" of FWM. In this way, the back scattering 
efficiency can be enhance. 

2.3.3 Phase Conjugate Beam Conversion 

For the output beam from an induction linac FEL amplifier or from a RF 
linac FEL oscillator without phase conjugate outcoupling, the transmission 
through the turbulent atmosphere would degrade the beam substantially. The 
FWM scheme described in the previous subsection can not be applied here 
because it requires two counterpropagating pump beam and the energy conver¬ 
sion efficiency is very low. For this purpose, a completely new nonlinear 
optics concept has been evaluated. It is a SBS process arranged in a FWM 
configuration (Figure 11) and, therefore, is called four wave stimulated 
Brillouin scattering (FWSBS). 

The system consists of a SBS cell placed in the path of the high power 
beam. A much lower power laser beam at the same frequency is incident from 
the opposite direction. The signal to be conjugated interacts at a very 
small angle with the two pump beams. The signal frequency is slightly 
higher than the pump frequency and the difference is equal to the acoustic 
wave frequency of the SBS medium. The result is a substantial energy flow 
from the high power beam into the direction opposite to the signal beam. 

The reflected wave is the phase conjugate of the signal field. This pro¬ 
cess has a high reflectivity (10 6 demonstrated) and high efficiency effect. 
Since the response time of the SBS process is in nanoseconds, it is appro¬ 
priate for the induction linac FEL. 

This new effect was studied in detail, both theoretically and 
experimentally, in this program. The results are discussed in Chapter III 
and IV respectively. 
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Figure 11 


Phase-conjugate beam conversion from a high power laser 
output, using the four-wave SBS effect 
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III. THEORY AND MODELING 


3.1 TASK OBJECTIVE 

The objective of the theory and modeling task in the FEL resonator 
development program is to develop analytical tools for the study of the 
nonlinear optics phenomena that have been identified for FEL applications 
and included in the conceptual resonator designs. To accomplish this 
objective, the approaches of formulating analytical theory and developing 
numerical modeling were followed for each phenomenon. The tasks performed 
are summarized 
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3.2 CONJUGATE FEL RESONATOR 

Conjugate resonators for general laser system have been an active 
research area for nearly a decade. Many analytical and numerical studies 
have been reported and were included in the reference list of two recent 
review articles 27 »28. However, most of the studies considered high Fresnel 
number cavities and assumed a large aperture conjugator with collimated 
input beams. In the conjugate FEL resonator design, the Fresnel number is 
low because of the highly focused beams. The fields at the conjugators and 
at the wiggler are in the far-field region of the field at the converging 
lens. The effects of aperturing on the conjugation process may not be the 
same for both resonator configurations. To facilitate the analysis of the 


01 - 127-86 


34 















effectiveness of conjugate FEL resonators, a two-dimensional steady-state 
simulation code was developed in this program. The code was used to anal¬ 
yze both conventional and conjugate FEL resonators. The results show the 
difference in mode structure and the improvement in aberration and align¬ 
ment control, when the conjugate reflector is used. 

3.2.1 FEL Resonator Modeling 

The phase conjugate FEL resonator code, FELPCR, was developed using 
wave optics and a simplified FEL gain model. The configurations include 
the conventional resonator with diverging lenses and the conjugator reson¬ 
ator with converging lenses. Schematics of both resonators investigated 
during the numerical studies are shown in Figure 12. 

The code uses cylindrical coordinates with radial grid point spaced 
for employing gaussian quadrature. The fields are propagated with propag¬ 
ators derived from the Fresnel-Kirchhoff integral. For convenience, the 
fields are decomposed into angular Fourier components. These components 
are coupled to each other at the tilted or aberrated mirrors. The aberr¬ 
ations are modeled by fifteenth-order Zernike polynomials. At the conjug¬ 
ator, the fields are complex conjugated assuming perfect conjugators. 

Since an RF linac FEL is essentially a low-gain device, the intracavity 
mode structures are dominated by the resonator configuration. The FEL inter¬ 
action provides simply the energy increase for the radiation, with an aper- 
turing effect due to the small electron beam cross-section. Based on these 
considerations, a simplified gain model has been developed for the FEL mode- 
medium interaction. It requires no particle simulation and generates a two- 
dimensional gain distribution on a single individual gain sheet. Considering 
the effect of betatron oscillation, the gain as function of radius is 
expressed as 29 


G(R) - 


/ g (Ur)) f (r)dr 


^ 0 f e (r)dr 


(2.1-D 
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Figure 12. Resonator configuration for resonator analys 
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where f e (r) is the electron density distribution, I(r) is the average 
optical intensity within the radius r, g(I) is the saturated FEL gain 
depending on the interaction intensity, I, 


g (i) 


'ss 


I + 


IJ3/4 

y 


( 2 . 1 - 2 ) 


In the above gain expression, g ss is the small signal gain and I s is the 
saturation intensity which is found to be near the resonance intensity for 
a tapered wiggler. This gain model is quite different from that for con¬ 
ventional lasers. There the gain is saturated by the local intensity and 
is thus, for example, lower in the high intensity region of the center of a 
gaussian optical beam. In the FEL gain model, however, the gain at a given 
radius is determined non-locally from the properties of all of the medium 
external to the given radius, properties which in turn depend on the aver¬ 
age inner intensity. The result is higher gain near the center of the 
beam. This peaking gain provides a natural aperturing effect for the FEL 
mode-medium interactions. 

3.2.2 FEL Resonator Analysis 

A comparison of the behaviors of both conventional and conjugate 
resonators was made using the code, FELPCR. The dimensions of the conven¬ 
tional resonator were chosen to be Li = 85 m, L 2 = 24 m and L e ff = 370 m. 
The effective cavity length, L e ff, is the distance between the end mirrors 
of the cavity which has the same mode size at the mirrors without using any 
diverging lenses. The choices of these values may not represent the actual 
design, but are largely limited by the Fresnel number that can be handled 
by the code. For the conjugate resonator, the converging lens is placed 30 
m from the wiggler and focuses the beam into the conjugator located 1.5 m 
from the lens. The beam size in the conjugator is approximately 50 fim. 

The mode structures are shown in Figure 13. The mode profile in a conju¬ 
gate resonator is close to a dome, instead of gaussian shape. This is 
consistent with the observations in other reports. 

The relative brightness has been obtained for different tilt angles 
and aberrations in both resonators. For the conventional resonator, the 
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Mode Intensity Profile 




Figure 13. Mode structuries in conventional and conjugate FEL resonators 



















misalignment sensitivity was found to be within 10 /<rad (Figure 14). In 
particular, the end mirror tilt is a factor of 5 more sensitive than the 
diverging mirror tilt. For the conjugate resonator, no obvious decrease in 
the brightness was observed up to 10 ^rad. Beyond this angle, the code 
showed oscillation in the output and can not converge properly, probably 
due to the small aperture size set at the conjugator. The clipper begins 
to interfere with the side-lobe field at the focus and causes the 
instability. 

For the study of the aberrations, we used the phase screen previously 
measured for NACL optics (Figure 15). The phase screen is scaled to obtain 
different rms values of the optical path difference (OPD). The effect of 
beam aberration was observed by applying the phase screen on the end mirror 
or the diverging lens of the conventional resonator, and on the converging 
lens of the conjugate resonator. It was found that the effect of aber¬ 
ration is more profound for the end mirror compared to the diverging lens 
(Figure 16). The aberrations have less effect on the performance of the 
resonator with phase conjugate mirrors. 

In general, we have found that the mode size in a conjugate resonator 
is larger than the mode size in a conventional resonator. The peak mirror 
loading is thus much lower. For alignment control tolerances, we observed 
that the sensitive tilt angle is well beyond 10 /<rad and might be near the 
mrad regime which is beyond the analytical capability of the code. The 
aberration control requirements were found to be relaxed by a factor of 
ten. The mirror surface requirement can be largely relaxed in a conjugator 
resonator. 

3.3 MULTIPLE-PULSE STIMULATED BRILLOUIN SCATTERING 

For backward SBS using a single short pump pulse, the reflected wave 
may not be observable because field does not have enough time to grow. For 
example, the pulse length from a RF linac FEL is typically 3 to 30 psec 
which is much shorter than the typical relaxation time of phonons (~ nsec) 
in most SBS media. However, if the next pump pulse arrives at the interac¬ 
tion region before the phonon field disappears, it permits the phonon field 
to grow (Figure 17). The result is the build-up of the phonon field from 
pulse to pulse. At steady state, the reflectivity near 100 percent even 
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RELATIVE BRIGHTNESS 



Figure 14. Cavity alignment sensitivity shown in brightness versus tilt 
angle for convertional (solid curves) and conjugate (dotted 
curve) FEL resonators. 
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Stimulated Brillouin Scattering 



phenomenon. The first observation was 
reflectivity of 96 percent. 























the Individual pulses are very short. This multiple-pulse SBS (MPSBS) 
phenomenon was first demonstrated at Garching using three consecutive 
pulses^l. The reported reflectivity was 96 percent. This effect can be 
applied to the phase conjugate FEL resonator and is studied in this 
program. 

3.3.1 MPSBS Theory 

The interaction of the pump, Stokes and phonons in a SBS medium can be 
described by the following one-dimensional equation,30 
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where E p , E s and Q are the slowly varying envelopes of the pump, Stokes and 
phonons, gj and g 2 are the coupling coefficients, v is the light velocity 
in the medium, T is the phonon relaxation rate. For most SBS media, the 
field Q grows slowly compared to the response time of the photon fields. 

In most cases, the optical fields are able to reach equilibrium within a 
few round trips of the medium length, which is of the order of 0.1 nsec. 
Therefore, the time derivatives in field equations can be neglected. In 
other words, for time scales longer than 0.1 nsec, it can be assumed that 
the fields response instantly to the change of phonon field, Q. 

An analytical solution of the Equation (3.1-1, 2, 3) has not been 
found with time dependence. However, if T is small enough it can be seen 
that the phonon field actually neglects the existence of the gap between 
the pulses. The SBS effect will be the same as if the pulses are packed 
together back-to-back. For the pulse length of 30 psec, it means that 
steady-stated can be reached in a couple of ten pulses. 

3.3.2 MPSBS Analysis 

In order to understand quantitatively the transient response of the 
SBS medium to short pump pulses, a one-dimensional code was developed based 
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on the equations (3.1-1, 2, 3). The code was designed to process every ten 
pulses as a group. At the beginning of each pulse, the Stokes intensity 
was assigned at noise level. Between pulses, the phonon field decays with 
the relaxation rate, T, (or relaxation time, r = 1/r). The reflectivity 
and the Stokes wave profiles are the output. 

Equations (3.1-1, 2, 3) can be converted into simpler forms in terms 
of a dimensionless quantity, G = gIL/r, where g = gig2i I is the pump 
intensity and L is the interaction length. A parametric study was per¬ 
formed to study the effect of relaxation time on MPSBS. The medium and the 
pump intensities were chosen such that G = 750/nsec. The pulse length was 
selected to be 0.1 nsec with a repetition rate of 2 GHz. For a gaseous 
medium, the relaxation time can be adjusted by changing the pressure. The 
interaction gain, g, is proportional to the relaxation time. The cal¬ 
culated variation of the reflectivity from pulse to pulse for different 
values of r are shown in Figure 18. The reflectivity was found to reach 
near the steady-state value within three pulses, which is consistent with 
the experimental observation reported earlier 2 *. The steady-state value of 
reflectivity, however, increases with the relaxation time primarily due to 
the increasing SBS gain. 

The reflected pulse profiles for the case of r = 1 nsec were examined 
and are shown in Figure 19. The pump intensity and length are normalized 
to one. In general, the pulse profiles have a sharp-rising front and fast- 
decreasing tail. The difference in steady-state reflectivity is mostly due 
to the difference in time for the pulse to reach the maximum intensity. 

For practical reasons, the case of a RF linac FEL using a xenon SBS 
cell was studied. For xenon gas at 39 atm, the SBS gain is 0.044 cm/MW and 
the relaxation time was 2 nsec 2 *. Assuming that the medium length is 1 cm, 
the pulse length is 30 psec with a repetition rate of 0.5 GHz. The case 
here shows that the relaxation time is equal to the pulse-separation time 
but several hundred times longer than the pulse length. The pump intensity 
varies from 10 to 100 GW/cm 2 and the noise level is chosen to be 0.1 
MW/cm 2 . The actual noise intensity may be much lower than this value. The 
change of intensity below the chosen level can be estimated by taking the 
spatial average of Eq. (3.1-2) assuming E p constant. For an intensity of 
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10 GW/cm 2 , the Stokes intensity doubles after every pulse. For high 
intensities, the rate should be proportionally increased. Thus, for every 
three orders of magnitude reduction in the noise level, the required number 
of pulses grows by ten. 

From Figure 20, It can be seen that the SBS reflectivity can reach its 
steady-state value within 10 or 20 pulses. For the cases of 50 GW/cm 2 , the 
reflected pulse profiles are shown in Figure 21. The first three reflected 
pulses are too weak to be seen in the linear scale. As the pulse number 
increases. The pulse profiles becomes closer to the shape of the input 
square pump. The steady-state reflectivity is 93 percent. 

In summary, the MPSBS effect has been studied numerically in real time 
domain. Choosing the RF linac FEL pulse structure and a practical SBS 
medium, it was shown that high reflectivity can be obtained after a small 
number of pulses. The reflectivity is close to unity that is high enough 
for FEL resonator applications. The reflected pulse shape is also close to 
the input profile. The FEL interaction is thus less perturbed. 

3.4 SHORT-PULSE FOUR WAVE MIXING 

The four wave mixing effect was considered as a phase conjugate 
outcoupling scheme for RF linac FEL resonators (Figure 22). For the anal¬ 
ysis of this phenomenon in the FEL characteristics regime, it was found no 
analytical theory and numerical modeling are available in the literature. 
Most of the theories presented the solutions either in steady-state or in 
long-pulse regime where the time varying pump fields are assumed to be 
uniform across the medium 32 . However, the pulse length of a RF linac FEL 
device is around 0.1 - 1 cm, which is probably shorter than the interaction 
region. The phenomenon is definitely in the transient regime and requires 
a full space-time representation in the coupled equations. In the area of 
numerical modeling, only the one-dimensional steady-state solutions have 
been reported. For the purpose of demonstrating the conjugation, field of 
view, beam misalignment and self-focusing effects, a two-dimensional simu¬ 
lation code is necessary to accomplish the task. The analytical theory and 
the numerical code developed in this program are described in the next two 
subsections followed by some detailed analysis. 
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Figure 20. Reflectivity build-up for the case of xenon SBS medium and 
typical short-pulse RF linac FEL at different intensities 
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Figure 22. FWM as an outcoupling scheme for FEL resonators 



























































3.4.1 SPFWM Theory 

The details of the analytical theory for SPFWM are shown in Appendix 
A, which is a preprint of the paper submitted for publication. The theory 
starts with the coupled-mode equations describing the four wave interac¬ 
tions. The equations can be solved exactly in the no-pump depletion 
regime. For arbitrary pulse profiles, the result is given in Equation (A- 
13). For a special case where the three input beams have the same pulse 
profiles, a closed-form reflectivity expression was obtained as shown in 
Equation (A-17). The result is quite different from the steady-state or 
long-pulse solution in earlier theories^, 32 # The reflectivity is indepen¬ 
dent of the pulse shape and monotonically increases with the pump fluences 
(Figure A-2 and A-3). 

Since bandwidth is one of the major issues for the FEL device, the 
theory is generalized to the cases where the coherence length is much shor¬ 
ter than the pulse length. The main results is in Equation (20), where 7 
is also the ratio of bandwidth to the transform-limited bandwidth. The 
dependence of reflectivity on this ratio is shown in Figure A-4. 

Consider a practical example for RF linac FEL pulses at 0.5 /tm 
wavelength. The pulse is 1 cm (i.e. 30 psec) and the pulse energy density 
is 10 J/cm 2 . If the bandwidth is one percent, the coherence length is 0.05 
mm. For a nonlinear optical medium similar to CS 2 , the reflectivity is 
found to be 

R = 200 [ I 0 2 ( 3 ) + Ij 2 (3) - 1 ] ~ 9000 ! (4.1-1) 

which is readily adequate for our requirement (> 10 3 ). 

The reflectivity expression was derived assuming coll inear polarization 
of the three beams. In this case, the equation for the conjugated field can 
be written as 


E 4 ~ *1111 E 1 E 2 E 3' (4.1-2) 

where prime means spatial derivative and ,nm is the nonlinear optical 
susceptibility. For circularly polarized input beams, as in our experi¬ 
ments, the x and y components of E 4 have contributions from different 
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combinations of the three input field components. Based on the nonzero 
third-order susceptibilities, the component E 4 X has the following 
contributions, 


E 4x ~ *1111 E lx E 2x E 3x + ^1122 E lx E 2y E 3y 

(4.1-3) 

+ *1212 E ly E 2x E 3y + *1221 E ly E 2y E 3x* 

For CS 2 , the off-diagonal suscepalibilities are related to ynn as 

*1122 ~ *1212 “ *1221 “ *1111 /<f (4.1-4) 
Since E ix = E 1y = E i for i = 1, 2, 3, we have 

E 4x " 3 4 ^ *1111 E 1 E 2 E 3 “ 1,1 *1111 E 1 E 2 E 3 f 4 - 1 " 5 ) 

Similar result should be obtained also for E 4 y. In the low reflectivity 
limit, the return intensity is thus about 2.4 times larger than the value 
in the col linear polarization case if the input intensities are the same in 
both cases. 

3.4.2 SPFWM Modeling 

For the purpose of studying numerically the SPFWM effect, a two- 
dimensional simulation code, called TFWM, has been developed in this pro¬ 
gram. 

In the four wave mixing configuration, beam 1 and 2 are the colinear 
pump beams propagating in opposite directions. Beam 3 is the probe beam 
incident at a small angle with respect to beam 1. The plane defined by 
these three beams is the two-dimensional interaction plane for the system 
and no variation of physical parameters is assumed in the direction normal 
to this plane. The z axis is defined as the line bisecting the angle bet¬ 
ween beam 1 and 3 (Figure 23), although the incident angle of each beam can 
be input independently. The return signal (beam 4) is expected to pro¬ 
pagate on the same plane and in nearly opposite direction of beam 3. 

The region of the nonlinear medium extends from 0 to L in the z 
direction and has an area large enough to cover the four beams. At the 
beginning (t = 0), the pulse fronts of beams 1, 2 and 3 are placed against 
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the front and rear edges of the medium. The advance of the pulses into the 
medium can be adjusted if a pulse delay is assumed in any one of the beams. 
For longitudinal sampling, the pulses are sliced into equally spaced sheets 
perpendicular to the z-axis. The slice distance should be twice that of 
the propagation distance in each time step, such that the field on each 
sheet will interact with every field sheet in the opposite-direction pulse. 

Transversely, the fields are sampled in equally spaced grid points 
across the whole aperture. The beam area is defined by either a gaussian 
profile or a smoothed square window function. The smooth window edges 
reduce the diffraction effect. Because the field sheets are not perpendi¬ 
cular to the actual direction of beam propagation for non-zero incidence 
angle, a linear phase variation is applied to the field across the beam. 

For short pulse applications, several different longitudinal pulse profiles 
have been used in the study (Figure 24). 

Before the pulses enter the interaction region, a phase screen 
representing the phase aberration can be applied to the complex field. The 
phase screen consists of five sinusoidally varying functions. With differ¬ 
ent variation amplitudes and periods, various phase screens can be created 
independently for each beam. The phase screen of beam 3 should be applied 
to the exiting beam 4 to check how well the aberration is corrected. In 
each case, the beam quality is calculated in terms of the Strehl ratio. 

The pulses advance in the medium using the near-field propagation 
scheme. In each propagation step, the one-dimensional field array is 
Fourier-transformed. For CW cases, the propagation continues until a 
steady state is reached. For pulse interaction, the propagation stops when 
the trailing edge of every pulse completely leaves the medium. At the end 
of the run, the reflectivity and Strehl ratio are calculated. The field 
arrays are stored for three-dimensional plots of pulse intensity and phase 
profiles. 

In the nonlinear medium, every field interacts with itself and with 
each other through the nonlinear optical susceptibility. The code models 
such interactions in Kerr-like media as well as in saturable absorbing (or 
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gain) materials. In a Kerr-like medium nonlinearities can be described by 
an effective index of refraction n: 

n = n 0 + n 2 <E 2 >, (4.2-1 

where <E 2 > is the time average of the square of the total field. Before 
each propagation, the complex field array, Aj (n x , n z ), are updated 
according to the formula 

Mj ■ -1 [ <2I 0 - IjlAj ♦ 2A*A m A n ] (4.2-2; 

J 

where the indices are 

j k m n 
12 3 4 

2134 (4.2-3) 

3 4 12 

4 3 12, 

and 

Ij * lAjl 2 , (4.2-4) 

I 0 = Ii + 12 + 13 + *4 
i = “ 

k 0 is the vacuum wave number, 0j is the incident angle of beam j and Az is 
the propagation distance. 

Equation (2) can also be transformed into coupled intensity and phase 
equations. Both methods are available in the code. Although using coupled 
intensity equations guarantees energy conservation, it is necessary to be 
very careful in the region where the intensity is zero because the inten¬ 
sity variable appears in the denominator of the accompanied phase 
equations. In general, the first term in Equation (2) represents the self 
phase modulation which would result in nonlinear self-focusing for a non- 
uniform beam. The second term represents the four wave interaction which 
generates the conjugated return signal. 




















3.4.3 SPFWM Analysis 

The two-dimensional code, TFWM, was developed for the observation of 
phenomena that have the nature of transverse field variations. The pheno¬ 
mena of interest related to SPFWM include self-focusing, effective field of 
view, conjugation fidelity and the saturation effect. 

In a nonlinear medium, self-focusing and FWM should have comparable 
effects because both are the results of electric field interactions through 
the third-order optical susceptibi1ity 3 ^. j n order to make sure that the 
self-focusing effect is properly included in the code, the process is 
isolated from FWM by inputing only one beam. The beam profile in the 
medium is carefully followed. At intensities above several GW/cm 2 , the 
self-focusing effects occur. A typical beam profile is shown in Figure 25 
for peak intensity of 2.5 GW/cm 2 . It is important to point out that the 
self-focusing observed is in only one transverse direction. It might be 
quantitatively different from the two-dimension self-focusing effect. For 
our purpose, one-dimensional self-focusing is sufficient and was well simu¬ 
lated by the code. 

Most analytical theories and one-dimensional codes of FWM considered 
only col linear interacting beams. However, the probe beam was incident at 
small angle with respect to the pump beam in most experiments and also in 
our applications. It is important to determine the dependence of the 
reflectivity on the incidence angle. The angular dependence was studied 
for the low reflectivity case so that it could be compared with the theore¬ 
tical prediction. The medium used was CS 2 in a 1 cm long cell. The length 
and the spatial width of the pulses were 1 cm and 0.25 mm, respectively. 

The pump and the probe intensities were 124 MW/cm 2 and 12.4 MW/cm 2 respec¬ 
tively. As the incidence angle was increased, the reflectivity decreased 
nearly linearly (Figure 26). The results were agreement with the dotted 
curve which was calculated considering that the reflectivity is 
proportional to the beam overlap. It was found that the beam-overlap 
argument was valid in most cases. 

FWM has been demonstrated experimentally as having a high conjugation 
fidelity. However, the FWM phase conjugation process has not previously 
been studied by a simulation code because no two-dimensional code was 
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Figure 25. Narrowing of pulse intensity profile (self-focusing) 
in a kerr medium (CS 2 )• The peak input intensity is 
2.5 GW/cm2 
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Figure 26. Angular dependence of reflecti 






















available. For demonstration, the intensity and the phase of an aberrated 
probe beam are shown in Figure 27. The aberration represents a Strehl 
ratio of 0.77. The beam interacts with two good quality pump beams in the 
CS 2 medium. The intensity and the phase of the reflected beam after pass¬ 
ing through the aberrator again are shown in Figure 28. The intensity 
profile is found to rise toward the trailing edge. This effect can be 
understood since the FWM energy transfer efficiency is proportional to the 
square root of the intensity product (liI 2 I 3 I 4 )• Since (I* + I 3 ) and 
(12 + I 4 ) are constants of the interaction, the intensity product will 
become higher as I 3 and I 4 increase. Therefore, the energy transfer is 

more efficient near the end of the pulse. The phase of the output field is 

uniform in most of the area. The correction of the aberration is nearly 
complete and the final Strehl ratio is calculated to be 0.997. 

The code is also used to demonstrate high reflectivities and the 
saturation effect. Figure 29 shows the theoretical prediction (solid 
curve) based on the constant pump approximation, and the simulation results 
(dots). It is demonstrated that high reflectivity can be obtained in the 
SPFWM effect. However, the output beam quality begin to degrade when satu¬ 
ration occurs. For 17 = 5, where the saturation is first observed, the 

output fields with square and gaussian input pulses are shown in Figure 30 
and 31, respectively. The intensity profiles are rich in structure and 
peak near the end of the pulse. The Strehl ratio of the reflected beam is 
found to be 0.84 and the energy transfer efficiency is 9.4 percent at a 
reflectivity of 800. 

3.5 FOUR-WAVE STIMULATED BRIL10UIN SCATTERING 

The concept of using four wave mixing in a RF linac FEL resonator was 
found to be useful in outcoupling the high power beam that is the phase 
conjugate of a suitable weak probe signal. However, it was found imprac¬ 
tical for FWM to be used for the induction linac FEL device due to the 
following reasons: First, although the FWM process requires two 
counterpropagating pump beams, there is only one high power output beam 
available. Second, the energy transfer efficiency in FWM is low, which is 
acceptable in a resonator because the remaining energy will continue to 
circulate. For the induction linac FEL, the non-transferred part of the 
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Conjugated Pulse After Aberrator 



nd phase of the conjugated and corrected 
aberrated probe beam shown in Figure 27 
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Figure 29. The saturation effect observed in simulation (dots) compared 
to the result from the non-depletion theory (solid curve) 
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beam energy would be lost. For successful induction linac FEL 
applications, a nonlinear process must have the properties of high 
reflectivity, high efficiency and require only a single high power beam. 

We have found a new nonlinear optics concept that meets these requirements. 
The new process, called FWSBS, was shown to be extremely efficient in 
converting the high power beam into the phase conjugate of a selected weak 
probe beam23,36. Both theoretical and experimental studies have been 
carried out on this program. 

3.5.1 FWSBS Theory 

In the FWSBS process, two pump beams, one at high power and one at low 
power, and having the same frequency, u, counterpropagate and meet in a SBS 
cell (Figure 32). The probe beam of frequency, u + Q, where Q is the 
phonon frequency, is incident at a small angle with respect to the high 
power pump beam. The arrangement is similar to the ordinary nondegenerate 
FWM configuration. An acoustic grating is created initially by the SBS 
interaction between the probe beam and the weak pump beam. This grating is 
further enhanced by the stimulated scattering of the high power pump beam 
in the opposite direction of the probe beam. The scattered field is a 
conjugate of the probe signal. 

No detailed theoretical study, especially concerning the saturation 
effect, can be found in the literature^?. since high reflectivity and high 
efficiency of FWSBS is of primary interest for the induction linac FEL 
applications, we have developed an analytical theory of FWSBS in this pro¬ 
gram. The detail of the theory is shown in Appendix B, which is a preprint 
of the paper submitted for publication. 

In the theory, the steady-state coupled equations of the four waves 
were formulated and solved analytically. The coupled equations (7) show 
that there are two kinds of interactions responsible for the FWSBS effect. 
The four wave mixing through the phonon interaction generates the initial 
return field. The mixing process ensures that the return field is the 
conjugate of the probe signal. As the return field becomes sufficiently 
strong the SBS amplifying process will take over and saturate the return 
field. The reflectivity and the efficiency in the nonsturation and 
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Four-Wave Stimulated Brillouin Scattering 



Figure 32. Schematic of four wave stimulated Brillouin scattering 
















saturation regimes were found to be given by the approximate expressions, 
which describe the dependence quite well, particularly when I 3 is less than 
12 . 


h + ! 3 


V ~ 1 + 1 ■ 
1 2 j 3 


(5.1-1) 

(5.1-2) 


where Ij, I 2 , I 3 are the intensities of strong pump, weak pump and probe 
beam, respectively. With properly chosen input intensities, both the 
reflectivity and the efficiency can be very high. These relations are 
compared to the experimental observations. No simulation code has been 
developed for FWSBS in this program. 

3.6 SYSTEM REQUIREMENTS 

Based on the analytical results obtained in previous sections. System 
analysis were performed for an up-link FEL system. The system is either an 
RF linac FEL resonator with FWM outcoupling or an induction linac FEL 
device with FWSBS beam conversion, the power requirements of the pump beam 
and the beacon can be readily determined from the required reflectivity and 
efficiency. In general, the results show that both nonlinear optics 
approaches are capable of accomplishing the atmospheric compensation for 
the high power FEL system. 

3.6.1 RF Linac FEL 

A four wave mixer in an RF linac FEL resonator can be used to 
outcouple the energy and convert it into the phase conjugate of a beacon 
signal which interacts with the intracavity pump beams. The fraction of 
the pump energy that goes into the conjugated wave is defined as the 
efficiency of the FWM process. From the resonator optics, it also repre¬ 
sents the outcoupling fraction. The interaction is assumed to take place 
in CS 2 or other media with comparable nonlinear optical coefficient. 

In the sort-pulse FWM phenomenon, the reflectivity was found to depend 
on the pump pulse fluence. It has a quadratic dependence in the low 
fluence region and increases exponentially at high fluences (See Figure 2 
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and 3 in Appendix A). Although it is highly underrepresenting at high 
reflectivities, the quadratic dependence is used in our analysis for sim¬ 
plicity. Figure 33 shows the relations of beam energies for different 
outcoupling fractions. For example, assuming 10 percent outcoupling frac¬ 
tion and 2 J/cm 2 pump fluence, the required probe fluence should be 20 
/*J/cm 2 and the reflectivity is 10 4 . 

Figure 34 shows the system requirement for different linewidths, 
assuming 10 percent efficiency. The lowest curve is for the case of trans¬ 
form limited pulses that are 3 psec long. The corresponding linewidth is 
0.1 percent at 1 pm wavelength. The curves move to higher pump fluences 
with increasing linewidths. Therefore, the effect of wider linewidth can 
be compensated by higher pump fluence in order to have the same efficiency. 
However, the reflectivity in this case is increased by the same factor. 

3.6.2 Induction Linac FEL 

The FWSBS interaction in an SBS cell can be used to convert the output 
beam from an induction linac FEL into a high power conjugate beam of a 
^ beacon signal. This is a high reflectivity and high efficiency process. 

If the power ratio of the two pump beams is sufficiently large, both 
reflectivity and efficiency are medium-independent quantities and are 
determined from the beam intensities only. Thus it is convenient to 
express the requirements in terms of intensity ratios (Figure 35). The 
relations were based on the approximate formulas in Equations (5.1-1) and 
(5.1-2). In general, the reflectivity can be very high depending on the 
pump intensity ratio. Efficiency near 50 percent should be easy to obtain 
if the probe intensity is close to the backward pump intensity. However, 
efficiency higher than 90 percent may be difficult because of the fast 
increase in the requirements of the probe intensity. 
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Figure 34. System requiremetns at different linewidths for an RF 
linac FEL using the FWM outcoupling scheme 
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Figure 35. System requiremetns for an induction linac FEL using the 
FWSBS beam conversion scheme 
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IV. EXPERIMENTS 


In the experimental portion of the program, experiments were performed 
to establish the viability of phase conjugators in FEL system configura¬ 
tions. Two four wave mixing phase conjugators were investigated. Liquid 
carbon disulfide was used to phase conjugate RF linac FEL-like short pulses 
to establish the scaling of reflectivity with pump intensity. Linear 
scaling of reflectivity up to 60% was observed for pump pulse fluence up to 
15 mJ/cm^. In the second series of experiments, SBS enhanced four wave 
mixing in acetone was used to phase conjugate induction linac FEL-like 
pulses to obtain reflectivities as high as 100 000% and 16% conversion 
efficiency. Both of these techniques show promise for the production of 
atmospheric aberration compensated high power beams in response to weak 
satellite beacons as described in Chapter II. 

This chapter consists of four sections. The first two describe how 
the two system concepts were studied experimentally, and how closely 
aspects of the experiments bear on those concepts. The two remaining sec¬ 
tions describe the four wave mixing (FWM) experiments in detail. The first 
section describes the short pulse four wave mixing experiments which used a 
doubled mode-locked Q-switched Nd:YAG laser source. Time averaged and time 
resolved reflectivity measurements are presented and analyzed. The second 
section describes the SBS enhanced four wave mixing (FWSBS) experiments, a 
type of non-degenerate four wave mixing, using Brillouin shifted pump 
beams. Measurements of reflectivity variation with probe, forward and 
backward pump beam intensities are presented. 

4.1 EXPERIMENTAL OBJECTIVES 

The objectives of the experiments were to evaluate the performance of 
high reflectivity/high efficiency phase conjugating mirrors used in FEL 
system concepts. Two types of four wave mixing phase conjugators were pro¬ 
posed in system concepts. In one concept, an extracavity phase conjugator 
is used to convert a high power pump beam supplied by an induction linac 
FEL into the phase conjugate of a weak beacon beam. High conversion effi¬ 
ciency is needed to avoid spoiling the overall system efficiency. In the 
second concept, an intracavity four wave phase conjugator is used to 
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outcouple an RF linac FEL, as well as provide compensation of atmospheric 
aberrations. Here, the conversion efficiency can be relaxed to the desired 
outcoupling, usually around 10% for high power RF linac FELs. Hence, an 
experimental objective was to measure efficiency. 

Our primary objective was to measure reflectivity scaling to anchor 
calculations of phase conjugate mirror reflectivities for high power FELs. 
In both system concepts, a weak beacon propagates through the atmosphere to 
the phase conjugating mirror, where the beam is time reversed with enormous 
gain. 

4.2 EXPERIMENTAL APPROACH 

The applicability of the experimental results derives from essential 
similarities between the FELs and our FEL simulators. For example, our 
experiments were performed in similar time regimes at similar wavelengths. 
However, the laser powers were much lower, and the linewidths much narrower 
than the FEL's. We sought to establish scaling relationships in power to 
predict performance at high powers. Also, our experiments did not facilit¬ 
ate investigation of the effect of large FEL bandwidth on the four wave 
mixers. Although attempts were made to expand the bandwidth of the mode- 
locked laser using self phase modulation, no significant broadening was 
achieved. Also, in attempts to emulate the induction linac FEL a second, 
relatively broad bandwidth laser source would not yield four wave mixing. 

(At the available power levels, narrow bandwidth pulses were necessary to 
achieve four wave mixing. Theoretically, increases in power can compensate 
for increase in bandwidth). Hence, the scope of the four wave mixing 
experiments was limited to the determination of reflectivity and efficiency 
scaling with pump and probe beam intensity variation. 

4.2.1 RF Linac FEL Phase Conjuqator Experiments 

To investigate the RF linac FEL intra-cavity phase conjugator, short 
pulse degenerate four wave mixing (DFWM) experiments were performed in a 
carbon disulfide; a fast Kerr-like medium. A mode-locked laser was used to 
produce a train of 25 ps pulses at 532 nm to simulate the optical pulse 
trains circulating inside the RF-linac FEL cavity. We chose to use CS? for 
the experiments because it is a Kerr-like medium a class which includes 
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other scalable, transparent, damage resistant candidates for high power 
devices, including liquids and plasmas. 

The short pulse DFWM experiments essentially replicate the geometry of 
the RF linac FEL intracavity four wave mixer. In the experiments, the 
three applied pulses are produced by a mode-locked laser. Note that in the 
system concept, the probe beam need not be a short pulse. A longer probe 
pulse would relax the difficult timing requirements set by the low duty 
cycle of the pump beams (~ 10"3). 

The FEL simulator for the short pulse (degenerate) four wave mixing 
experiments was a Q-switched, mode-locked doubled Nd:YAG laser. Table 7 
compares the intracavity FEL pulses to the FEL simulator output. Figure 36 
shows a schematic of the simulator output defining several time constants 
listed in the table. The laser emits a burst of pulses every 100 ms. 

Each burst consists of approximately 12 pulses with energies distributed 
under a Gaussian envelope. The pulses have an average duration, of 25 ps. 
(The uniformity of the pulse duration over the burst is discussed in 
section 4.3.3). 

Table 7 shows that the performance of the mode-locked laser is similar 
to an RF linac FEL in many respects, including <i>he pulse peak power, energy 
per pulse, pulse width, pulse spacing and wavelength. The pulse energy of 
5 mJ compares favorably with the minimum FEL pulse energy of 10 mJ; 
however, it must be borne in mind that the FEL pulse energy inside the 
cavity is about ten times as great, and that the pump pulse supplied to the 
four wave mixer is reduced by a factor of tour by losses in the optical 
train and division of the pulse into two pump pulse. Also, note that 
because a fast medium is used in the experiment, with a picosecond response 
time, the temporal profile of the pulse train is not important. The medium 
responds to each pulsed interaction as a separate event. 

4.2.2 Induction Linac FEL Phase Conjuqator Experiments 

To indicate the potential for the application of a four wave stimu¬ 
lated Brillouin scattering (FWSBS) phase conjugator to an induction linac 
FEL, Table 8 compares the nominal characteristics of the (induction linac) 
FEL simulator with the FEL. The "demonstration requirements" use 
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Table 7. PEL System Performance and Current Simulator Capabilities 


Parameter 

FEL 

Quantel 501 

power 

1-10 MW 

620 kW [1,3] 

peak power 

1 GW 

210 MW [1] 

energy/pulse 

10-100 mJ 

5.2 mJ [1] 

pulse width 

10-100 ps 

25 ps [2] 

pulse spacing 

2-10 ns 

7.1 ns [1] 

wavelength 

Visible 

532 nm [l] 

bandwidth 

0.5 - 1.0% 

<0.04% [1] 

Notes 

1. Measured 

2. Manufacturer's specification 

3. Averaged over the 6 most energetic pulses 


Table 8. Nominal Laser Characteristics for FWSBS 


Demonstration 

Requirement 

Does the Induction linac 

FEL Capability meet the 
demonstration requirement? 

Pulse Duration, > 15 ns 

Yes 

Peak Power > 2 MW 

Yes 

Beam Quality < 1.04 

Yes 

Wavelength = 1.06 fi m 

Yes 

Relative Bandwidth, hvlv < 2 x 10' 7 

No 




























Figure 36. Temporal structure of mode-locked, Q-switched Nd:YAG laser 
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the FEL simulator characteristics as the minimum requirements for FWSBS. 
For this type of four wave mixing, the FEL characteristics equal or exceed 
those of the simulator, with the exception of linewidth. As was the case 
for the DFWM experiments, reflectivity scaling with power may compensate 
for the much larger linewidth. 

As mentioned previously, FWSBS is a nondegenerate type of four wave 
mixing that requires a small frequency difference between the pump and 
probe pulses. In our experiment, we obtained this shift by using a Brillo- 
uin scattered pulse for the pump. However, an actual FEL implementation of 
a FWSBS phase conjugator would not pre-shift the pump beam wavelength as we 
did in our experiment. Instead, the small wavelength difference between 
the pump beams and the probe beam necessary for SBS enhancement can be 
obtained by design of the FEL and the beacon. This avoids the losses of an 
additional SBS cell in the FEL pump beam. 

The FWSBS investigation consisted of several reflectivity measurements 
and an experiment which optimized the pump beam conversion efficiency. The 
variations are shown in the test matrix presented in Table 9. For each 
matrix entry, reflectivity was measured as the probe beam intensity was 
varied over four orders of magnitude. This data was used to study the 
dependence of the reflectivity on 1 3 , Ij and 1 2 /1 1 where I], I? and I 3 are 
the forward pump, backward pump and probe beam intensities. In addition to 
these experiments, an experiment in which the pumped volume was enclosed by 
the probe beam was performed to maximize the pump energy extraction. 


Table 9. FWSBS Test Matrix 


I 2 /I 1 

I] 

4% 

1% 

0.048% 

46 MW/cm 2 

X 



80 MW/cm 2 

X 

X 

X 
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4.3 SHORT PULSE DFWM EXPERIMENTS 


This section describes the degenerate four wave mixing experiments 
using green picosecond pulses in carbon disulfide that were performed to 
study phase conjugation reflectivity scaling with pump pulse fluence. The 
experiments showed a power law dependence of reflectivity on pump pulse 
fluence, with reflectivities ranging up to 60% for pump pulse fluences up 
to 15 mJ/cm^. Fluence was predicted to be the correct scaling parameter, 
rather than intensity or energy, for short pulse four wave mixing in 
Appendix A. 

A power law dependence was observed in both pulse train average and 
single pulse experiments. An active-passive mode-locked Nd:YAG laser 
oscillator and amplifier were used to produce a pulse train of about 12 
pulses under a Gaussian envelope. Each pulse was split into pump and 
signal pulses and used in the four wave mixing geometry to obtain a phase 
conjugate reflection. In pulse average reflectivity experiments, the inte¬ 
grated pulse train reflectivity was measured with respect to the pulse 
train fluence. In the single pulse measurements, reflectivities were 
obtained from fast photodiode traces of the probe and phase conjugate 
reflection, with scaling by pulse train energy measurements. 

4.3.1 Experimental Set-Up 

The experimental set-up is shown in Figure 37. Bursts of short pulses 
are emitted at 10 Hz by the double mode-locked, Q-switched Nd:YAG laser (to 
be described in detail below). The pulses first pass through a polarizer 
and quarter wave plate which provide isolation to the laser. A 4% 
beamsplitter samples the pulses to form the probe beam while the pump 
pulse passes on to a spatial filter. This round-about path was chosen to 
allow the introduction of a self-phase modulator, (SPM), to increase the 
pump pulse bandwidth. Then, at BS, the pulse is split into two pulses 
which are made to cross paths at the four wave mixing cell (FWM). After 
the 4% beamsplitter, the probe beam is demagnified by a long focal length 
telescope, passed through a variable delay line, through a 50 ; ; beam sampler 
and on to an insertion mirror which steers the probe towards a “collision" 
with the pump pulses. 








































A mode-locked laser system, with peak powers as high as 1 GW, was used 
for these experiments. The system consists of a mode-locked Q-switched 
Nd:YAG oscillator and amplifier, followed by a second harmonic generator. 

As shown in Figure 36, the system emits trains of 10-12 short pulses (25 ps 
duration) at 10 Hz. The unit was modified by the removal of the standard 
pulse selector, so that the entire pulse train is output. The performance 
specification of our modified Quantel YG-501-10 are given in Table 10. 

Table 10. Quantel YG 501-10 performance specifications 


Parameter 

Performance 

Notes 

repetition rate 

10 Hz 


wavelength 

532 nm 


beam diameter 

7 mm 


polarization 

vertical 


pulse train energy, E 

2 mJ< E <25 mJ 

Variable with amplifier 
flashlamp voltage 

peak single pulse energy, 
Ed 

12% 

Photodiode trace analysis 

P 

pulses per mode locked 
trai n 

10-12 


period between pulses 

7.1 ns 


pulse duration 

25 ps 

Manufactures specification 

TPA measurement 

pulse train stability 

5% 


bandwidth 

<2.5 A 

Resolution 1imit of 
monochrometer 


>0.17 A 

Transform limited linewidth 


The pulse train energy is the total of the individual pulse energies 
in a train. The distribution of energy among the pulses was monitored by a 
photodiode. A typical photodiode signal is shown in Figure 40. The 
apparent fractional nanosecond width of the signal peaks is due to instru¬ 
ment broadening by the detector in combination with oscilloscope. In this 
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trace* the pulse energies may be inferred from the relative heights of the 
peaks. Typically, the peak pulse contains approximately 12% of the pulse 
train energy. 

No rigorous measurements of the pulse duration or bandwidth were made. 
The manufacturer claims a pulse duration of 25 ps, which we corroborated 
qualitatively with a two photon fluorescence measurement of the fundamental 
1.06 /im pulse width [40]. Limits on the bandwidth were established by the 
transform limit of 0.17A, determined from the pulse duration, and instru¬ 
ment limited linewidth measurements of 2.5 A. Also, a Michelson 
interferometer gave good fringe contrast over about 1 cm of path length 
difference, from which we estimate a bandwidth of about 0.27 A. 

An optical schematic of the laser system is shown in Figure 38 The 
oscillator is an active/passive mode-locked Nd:YAG laser [38,39]. 
Simultaneous mode-locking and Q-switching is caused by a cell containing a 
saturable absorber in the cavity. For this purpose, a flowing dye cell in 
contact with the rear mirror is used to circulate a solution of Kodak 9740 
Q-switching dye in 1, 2 dichloroethane. In the presence of sufficient dye 
concentration, the laser mode-locks and Q-switches, resulting in a burst of 
pulses spaced by the cavity round trip time of 7.1 ns. With increasing dye 
concentration, the energy of the burst increases, fewer pulses occur in the 
burst and the pulse width is shortened, with the net effect of increasing 
the second harmonic generation. The mode-locking properties of the dye 
tend to decay over a few days, necessitating regular dye changes for ener¬ 
getic pulses. 

The introduction of the acousto-optic mode-locker into the cavity has 
the effect of stabilizing the mode-locking. The pulse train energy stabil¬ 
ity was found to be approximately 5%. 

The pulses were mixed in an 8 mm thick spectrophotometer cell contain¬ 
ing spectroscopic grade carbon disulfide. The pump beams were made to 
focus at the cell by slightly defocussing the spatial filter, giving an 
elliptical spot of 1.05 mm x 1.31 mm. The probe beam was demagnified by 
about 3 and brought to focus at the cell using a 1 m lens (not shown), 
giving a spot size of 0.46 mm at the cell. The probe beam was introduced 
into the pump beam path at an angle of about 10 mr. 
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tel Mode-Locked Nd:YAG Laser System. 1: flowing 
; 2: acousto-optic mode-locker, 3: variable dela 
; 4: Nd:YAG amplifier; 5: oscillator output coup 
('AG amplifier; 7: second harmonic generator. 





























The coincidence of the three incident pulses was obtained interfero- 
metrically. First, to determine the location at which the pump pulses 
crossed, a front surface mirror was inserted at the approximate location of 
the crossing such that the two pulses were retroreflected. Then, using a 
beam sampler placed upbeam of the pump pulse splitting, the retro-reflected 
pump beams were projected onto a screen. When the front surface of the 
mirror located the pump crossing, interference fringes could be observed on 
the screen. Similarly, when properly timed, the inside edges of the for¬ 
ward pump and probe pulse could be made to interfere when the two beams 
were made nearly parallel. The timing of the probe pulse was adjusted 
using the right angle prism optical delay line. 

Data was collected at the three positions occupied by PD, D1 and D2. 
The output from PD, a fast photodiode, was used to display the mode-locked 
pulse train on a 500 MHz bandwidth osci1loscrope (Tektronix 7844), as in 
Figure 40. Because oscilloscope traces at 10 ns/division are faint, 
visible film records of pulse trains required as many as 200 multiple sweep 
exposures. Detectors D1 and D2 were pyroelectric detectors (accessories to 
the Laser Precision Rj 7200 energy meter) used to measure probe and signal 
pulse train energies. To perform time resolved reflectivity measurements, 

D1 and D2 were replaced by fast photodiodes. 

In order to measure the pulse energies, transmission factors with 
respect to the energy detector (the instrumentation is described below) 01 
were measured. Since the probe and pump beams were both proportion to the 
laser output, a measurement at D1 gave all of the incident beam energies. 

In addition, because the probe beam energy was proportional to the pump 
beam energy, the efficiency of the four wave mixing process was propor¬ 
tional to the reflectivity. 

In an effort to study the dependence of reflectivity on bandwidth, we 
attempted to modify the pump pulse bandwidth using self phase modulation 
[41]. As shown in the figure, the layout was designed to accommodate the 
insertion of a 20 cm cell of CSj. Self phase modulation was evidenced by 
several experimental results, but was difficult to quantify as discussed 
below. 
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4.3.2 Time Averaged Reflectivity 


Figure 39 shows a plot of probe pulse train reflectivity measurements 
versus pump pulse train energy. In the data, the pump pulse train energy 
refers to the sum of the pulse energies in one of the pump beams, = E 
Eii, where the Eii are individual pump pulse energies. The summation is 
performed by the pyroelectric detectors which cannot resolve the individual 
pulses in the pulse train. Similarly, the average reflectivity is given by 
R = (£ E 4 i)/ (E E 3 -j) where the £ 4 } and E 3 -j are individual signal and probe 
pulse energies, respectively. In the plot, the data is seen to have a 
powerlaw dependence on the pulse train energy, given by the slope of the 
line fit to the data, R a E 1*3. Another notable feature of the data is 
the sharp drop in reflectivity at fluences in excess of about 200 mJ/cm 2 . 
The drop may be caused by the onset of self-focusing in the CS 2 , or the 
onset of self phase modulation in the Nd:YAG amplifier rod with high pulse 
intensities, broadening the pulse bandwidth. 

For comparison, a plot of the theoretical reflectivity for negligible 
pump depletion is shown. The calculated pulse train reflectivity is cal¬ 
culated from R, = 1.8 rjj?, where rj] = 63*F-j (J/cm 2 ) (specific to CS 2 ) and 
Fj is the individual pump pulse fluence in units of J/cm 2 . (See the 
following section for discussion of this expression.) The disagreement 
between the slope of the curves is probably due to the significant pump 
depletion at these probe intensities. Although the probe pulse energy was 
less than 5 % of the pump pulse energy, the relatively small probe spot size 
gave the probe 27% of the pump pulse fluence. Certainly the peak pump 
pulses, with reflectivities in excess 65%, experienced significant (> 17%) 
pump depletion. 

4.3.3 Time Resolved Reflectivity 

Time resolved reflectivity measurements were used to obtain single 
pulse reflectivity dependence on pump intensity. Typical time resolved 
reflectivity data is shown in Figure 40. Traces of the probe and phase 
conjugate signals were recorded simu!taneously using photodiodes and a dual 
beam oscilloscope. The upper and lower traces correspond to the phase 
conjugate and probe signals, respectively. Notice that the envelope of the 
phase conjugate reflection rises and falls more sharply than the probe he 
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Conjugate return beam 


probe beam 



Figure 40. Typical time resolved reflectivity measurement. 

The upper trace corresponds to the phase conjugate 
- reflection, the lower to the probe pulse train. 

(20 ns/division). 


envelope, indicating a nonlinear reflectivity for pulses at the leading and 
trailing edges of the pulse train. 


Figure 41 Is a plot of individual pulse reflectivities as a function 
of pump pulse fluence. The numbers next to the data points refer to the 
pulse position in the pulse t.ain, showing a rise to a peak reflectivity of 
65% for the sixth pulse. The points are well fit by a straight line in the 
log-log plot, which has the nearly linear form R a The four wave 

mixing efficiency for each reflection as given by the right hand vertical 
axis of the figure. Since the probe was proportional to the pump beam is 
these experiments, the efficiency, or fraction of pump energy transferred 
to the phase conjugate beam, was proportionate to the reflectivity. To 
take the small size of the probe beam into account, only the portion of the 
pump beam overlapping the probe beam was considered. 


For comparison, the predicted reflectivity of a transform limited 
pulse in the non-depleted pump regime is shown. The theoretical curve is 
derived from the expression for the reflectivity given by eq. 17 of 
Appendix A. That expression considered the interaction of linearly 
polarized beams; the reflectivity must be increase by a factor of ~2 to 
discribe circularly polarized beams (as was the case in our experiment), as 
explained in Seciton 3.4.1. Hence, the dependence can be expressed as, R = 
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1.8 rj2, where 7 Is a parameter proportional to the nonlinear optical sus¬ 
ceptibility and the pump fluence. Using the jO) value of 2.025 x 10'31 In 
SI units for CS 2 , 7 is about 63 times the pump fluence in units of (J/cm 2 ), 
as given in Equation (19) of Appendix A. 

Although the experimentally measured reflectivity agrees with the 
theory for a pulse fluence of about 12 mJ/cm 2 , the data does not show the 
predicted quadratic dependence of reflectivity on pump fluence. Two plau¬ 
sible explanations for the discrepancy are: (1) Pump beam depletion 
reduces the reflectivity or ( 2 ) the pulse bandwidth was proportional to 
pulse intensity, reducing the reflectivity to a linear dependence on 
fluence. With a peak conversion efficiency of 27%, pump depletion cer¬ 
tainly affected the results. With respect to the second explanation, band¬ 
width modulation over pulse train produced by mode-locked Nd lasers has 
long been observed in passively mode-locked neodynium lasers [42,43] and 
recently measured in doubled active/passive mode-locked Nd:YAG laser system 
quite similar to ours by Nathel [44]. Nathel’s simultaneous time resolved 
measurements of pulsewidth and bandwidth showed that both quantities vary 
significantly over the duration of the pulsetrain. Up to the peak pulse, 
the seventh in a train of 13 pulses, the pulse train exhibited a nearly 
constant pulse duration (about 30 ps), but has bandwidths that Increased in 
proportion t the pulse energy. By the peak pulse, a four-fold increase in 
bandwidth over the first pulse has occured. After the peak pulse, the 
bandwidth stabilized with the accompaniment of a rapid decrease in the 
pulsewidth to just a few picoseconds. 

Such evolving pulse characteristics may explain our results. 

According to predictions of the theory (see Appendix A), the reflectivity 
should be inversely proportional to the bandwidth. With bandwidth propor¬ 
tional to the pulse energy and therefore fluence, the pulsetrain reflec¬ 
tivities are expected to be linear in pulse fluence. Pulse reflectivities 
after the peak may be degraded as the error in the pulse timing becomes 
increasingly significant with respect to the pulsewidth. 


4.3.4 Bandwidth Effects on Reflectivit 
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resolution capabilities. However, two indications of our success in' 
achieving nominal bandwidth expansion were observable effects on the time 
resolved reflectivity measurements and a narrowing of the four wave mixing 
reflectivity versus probe delay characteristic. 

We attempted to use self phase modulation (SPM) to control pulse 
bandwidth. Self phase modulation, (SPM), the process responsible for the 
dynamic bandwidth modulation in mode-locked neodymium lasers, occurs when 
an intense transient pulse propagates through a medium with non-linear 
refractive index 03 , such as Nd:YAG (n 2 = 4 x 10 -13 esu) or CS 2 (n 2 - 1 x 
10“1 1 esu) [45]. The dependence of the index of refraction on the pulse 
Intensity causes portions of the pulse to propagate with differing phase 
velocities. For 03 > 0 (as in CS 2 ), this has the effect of red-shifting 
the leading edge of the pulse and blue-shifting the trailing edge, giving 
the pulse a frequency sweep or “chirp." We performed experiments with a 20 
cm cell of CS 2 introduced into the path of the pump beam to investigate the 
properties of DFWM with an SPM-broadened pump pulse and relatively 
narrowband probe pulse. However, because we could not broaden the pump 
pulses beyond the resolution limits of our monochrometer (2.5A), we were 
not able to quantify bandwidth dependence. 

The effects of the SPM are demonstrated by the time resolved 
reflectivity oscillograms shown in Figure 42. The figure compares the 
phase conjugate signal (upper trace) for pump pulses in and out of the 
regime of s 5 gnificant self phase modulation. In (a) and (b) the intensity 
of the peak pulse at the SPM was about 310 MW/cm? and 470 MW/cm?, respec¬ 
tively. In (b), SPM broadening of the peak pump pulses causes the suppres¬ 
sion of the phase conjugate signal. 

The broadening of the pump pulse bandwidth by the SPM is also evident 
is measurements of the reflectivity with probe delay. Leading and trailing 
edges of the pump beam which have been frequency shifted by the SPM destroy 
the degeneracy of the four wave mixing and spoil the reflectivity. Figure 
43 shows measurements of the reflectivity with probe pulse delay, with and 
without the self phase modulated pump pulses. The characteristics have 
widths of 25 and 19 ps, for measurements with and without the self phase 
modulator, respectively. 
















Conjugate return beam 


probe beam 



Conjugate return beam 


probe beam 


Figure 42. Effect of self phase modulation on time resolved reflectivity 
measurements. The two oscillograms show time resolved 
reflectivity measurement with the self phase modulator in the 
pump beam path. In (a), the peak pulse intensity at the SPM 
reaches about 310 MW/cm 2 and in (b) 470 MW/cm 2 . At the 
higher pump pulse intensities (b), the broadening becomes 
significant, and the reflectivity is reduced. 
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Figure 43. Reflectivity versus probe pulse delay with and 
without the self phase modulator. 
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4.3.5 Conclusion 

In the short pulse DFWM experiments, we were able to observe phase 
conjugate reflectivities as high as 65%, the highest value ever observed 
for short pulse DFWM as far as we know (compare Reference 48), and an 
efficiency as high as 27%. The observed scaling of reflectivity with pump 
pulse fluence did not obey the quadratic scaling law predicted by the 
theory, but this discrepancy is most likely due to variations in the pulse 
bandwidth in the course of our experiments. 

The difficulty imposed by the bandwidth variation over the pulse train 
can be remedied by the use of a single pulse selector. By selecting the 
same pulse from each pulsetrain, fixed bandwidth experiments could be per¬ 
formed. 



Bandwidth control might be achievable using schemes other than self¬ 
phase modulation. One possibility would be to broaden single pulses by 
propagation through long single mode fibers followed by a relatively broad¬ 
band Nd:g1ass amplifier. The use of optical fibers to introduce chirp into 
short pulses Is part of a widely used technique for pulse compression. 
Another approach might be to use the doubled output of the mode-locked 
laser system to synchronously pump a dye laser with intra-cavity bandwidth 
tuning elements. 


4.4 FOUR WAVE STIMULATED BRIILOUIN SCATTERING 


This section describes and analyzes the Four Wave Stimulated Brillouin 
Scattering (FWSBS) experiments which produced reflectivities as high as 
100,000'c and efficiencies as high as 16 c ;. First, the experimental set-up 
and its characteristics are described, followed by presentations of reflec¬ 
tivity and efficiency measurements. Reflectivity measurements with varying 
forward pump, backward pump and probe intensities are presented and 
compared to simple theoretical predictions. The results of an extraction 
efficiency experiment are presented and compared to a simple theoretical 
prediction. 


4.4.1 Experimental Set-Up 

The experiment was arranged as shown in Figure 44. A Q-switched 
Nd:YAG laser system (described in detail below) produced a 12 ns pulse 
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which was first sampled at BS to produce the probe beam and then back-scat¬ 
tered by the acetone SBS cell (SBSC) to form the strong forward pump beam. 

The pump beam polarization was turned 90° by two passes through a quarter 
wave plate and reflected out of the retrobeam path by the polarizing beam¬ 
splitter (PBS). The strong pump is then focused by a 3 m focal length 
meniscus lens to an elliptical spot of 1.5 mm x 2.3 mm at the acetone four 
wave mixing cell, FWMC. After passage through the cell, the beam encoun¬ 
tered the normally oriented surface of an uncoated wedged beam splitter 
which provided the weak backward pump. After being split off the main 
beam, the probe beam was sampled by a 50% beamsplitter, focused by a 1 m 
focal length lens and clipped in to the strong pump beam path at an angle 
of about 10 mr. 

We used a Q-switched Nd:YAG laser system operating at 10 Hz to produce 
pulses in a horizontally polarized beam, which we clipped to a roughly 
circular cross-section 7 mm in diameter. The laser system, which consisted 
of an oscillator and amplifier, was capable of producing pulses of up to 
250 mj, although we were limited to energies of 150 mJ by the onset of 
damage to our Brillouin shifter. A typical temporal profile of the pulse 
as shown in Figure 46. The pulses were smooth with widths of approximately 
12 ns as long as single axial mode operation was maintained by the oscil¬ 
lator. Well tuned single axial mode operation could be maintained for 
minutes at a time, before modulation in the pulse envelope had increased 
noticeably. Hence, the laser was frequently fine tuned to minimize modu¬ 
lations in the pulse envelope due to the presence of additional mode com¬ 
ponents. This instability contributed to scatter in the experimental data. 

The laser pulse was Brillouin shifted upon backscattering by the SBS | 

cell (SBSC) to obtain the pump pulse at the desired wavelength. The SBS 
cell consisted of a 8" x 1" l tube with optical windows expoxied to its 
ends. The cell was half filled with spectrophotometric grade acetone, and 
topped with argon gas to avoid combustion of the acetone vapor in air. The 
cell is shown in the plane of the optical bench for purposes of the 
schematic; in fact, the cell was oriented with its axis slightly off the 
vertical. This allowed the f,20 beam to enter the cell at a relatively low 
intensity, avoiding damage to the cell window, while reaching focus at 
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S' about 2 cm below the surface of the acetone, which optimized the SBS 

reflectivity. Typically the SBS reflectivity was about 70%. 

The four wave mixing cell, labelled FWMC in Figure 44, consisted of a 
4 cm x 1" ^ cell, capped by optical windows and filled witl spectrophotc- 
metric grade acetone. The cell axis was always misaligned with respect to 
the pump and probe beams to avoid multiple weak backward pump beams and to 
avoid retroreflecting the probe beam. 

Energy and pulse profile measurements were made by pyroelectric probes 
(accessories to the Laser Precision Rj 7200 energy meter) Dl, D2, D3, D4 
and photodiodes PD1 and PD2 (PD3 was used to trigger oscilloscopes). Dl 
and D2 were used to monitor the SBS and FWSBS reflectivities, respectively. 
The photodiodes were used in conjunction with Tektronix 7104 oscilloscopes, 
which were capable of producing bright enough traces to record nanosecond 
signals on film. PD2 was used to monitor the shape of the probe and phase 
conjugate signals. PD1 was used to monitor the shape of the forward pump 
beam before and after transmission through the four wave mixing cell. 

Since the backward pump was produced by sampling the forward pump pulse 
^ after transmission by FWMC, the second pulse signal monitored by PD1 may be 

viewed as either the backward pump pulse or the post four wave mixing for¬ 
ward pulse. 

Because thermal lensing of the forward pump pulse into the phase 
conjugate signal path occurred for operation at 10 Hz, a pulse selector was 
used, shown as PS in Figure 44. At 10 Hz operation, a false "phase con¬ 
jugate" signal would grow over the course of a few seconds with the probe 
beam blocked. This occurred as the pump beam heated the the acetone, and 
thermally induced refraction steered scattered rays into the signal beam 
path. To avoid this, we used an electro-mechanical shutter, synchronized 
with the laser flashlamp timing to pass a single pulse upon command. By 
waiting about 10 seconds between shots, the false signal with the probe 
beam blocked was eliminated. 

4.4.2 FWSBS Reflectivity Experiments 

FWSBS reflectivity measurements were made while varying pump and probe 
intensities. Figure 45 summarizes data collected from four reflectivity 
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experiment. Each point represents a single four wave mixing event. In 
each of the plots, the reflectivity variation with probe intensity is 
shown. Note that is comparing (a) and (b), only the pump power was 
changed. Plots (b), (c) and (d) differ only in the magnitude of the 
backward pump beam relative to the forward pump beam. 


Each of plots includes a theoretical reflectivity given by: 


! 2 + ! 3 


(B-27) 


where Ij, I 2 and I 3 are the intensities of the forward pump, backward pump 
and probe beam, respectivity. (See Appendix B for the derivation of Equa¬ 
tion B-27). This characteristic asymptotically approaches the value of 
I 1 /12 for small values of probe intensity, and falls as I/I 3 for I 3 > I 2 . 
The theory appears to characterize the data well in the regime where I 3 J, 
I 2 ; however, much larger reflectivities than were predicted by the theory 
were obtained for small I 3 . 

This discrepancy between the theory and the data can not be explained 
in terms of signal noise. Background noise was measured for the data in 
plots (a), (b) and (d), by blocking the probe and measuring the signal in 
the return beam path. The noise level during the I 2 /I 1 = 4% experiments 
was stable and small, and has been subtracted out of the phase conjugate 
signal energies. In plot (d), where the backward pump intensity was 
reduced by two orders of magnitude, the background noise was significant 
for probe intensities under 1 MW/cm^. However, because the scatter in the 
background noise measurements was large, no systematic error could be 
discerned. In plot (d), the four smallest reflectivities for I 3 = 0.5 
MW/cm^ and 0.05 MW/cm^ may be entirely due to noise. 

Typical FWS 8 S pulse profiles are shown in Figure 46. It is seen that 
the return pulse is about 50% narrower than the probe pulse. Because the 
FWSBS medium is short compared to the pulse lengths (in space), the shape 
of the return pulse envelope is related to the product of the envelopes of 
the three incident pulses. Note that the backward pump pulse, which is 
also a sample of the forward pump pulse after transmission through the 
acetone cell, has been considerably shortened. This shortening may 
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Figure 46. Typical FWSBS pulse profiles. 
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correspond to significant depletion of the trailing edge of the pulse after 
the onset of SBS enhancement of the four wave mixing. The risetime of the 
SBS enhancement is approximately the phonon lifetime in acetone of 2 ns. 

A summary of the experimental Investigation into the dependence of 
reflectivity on pump power is shown in Figure 47. In the figure, the data 
presented in Figure 48(a) and (b) is replaced by averages and error bars, 
and fit using linear regression. The error bars are constructed such that 
they enclose all of the data points at a fixed probe beam intensity. The 
reflectivity roughly increased by a factor of four for a two-fold increase 
in pump power. This scaling is faster than the rate given by equation B- 
27, which predicts a proportional increase in R. 

A summary of the experimental investigation into the dependence of 
reflectivity on the relative intensity of the backward pump beam is given 
in Figure 48. The plots summarize the data shown in Figure 45 (b), (c) and 
(d) in a manner similar to the previous figure. The line fits exclude the 
weakest probe intensity data points for each curve, which are outside the 
pump depletion regime where R a I/I 3 . 

From the plots, it may be observed that the reflectivity increases 
with I 2 /I 1 . contrary to the trend predicted by equation B-27. This 
disagreement with the theory is to be expected in the limit of small I 3 . 
While the theory predicts a diverging reflectivity in the absense of the 
backward pump as I 2 ■* 0 , the elimination of the backward pump would 
extinguish the interaction. 

The diminution of the reflectivity with the backward pump may be 
understood as a transition from a well defined backward pump to a diffuse 
backward pump. As the backward pump intensity approaches the intensity 
level of randomly scattered pump beam in the cell, the scattered light 
could begin to play the role of backward pump. This noise in the backward 
pump would be imprinted by the four wave mixing on the outgoing signal. 

The fidelity of the phase conjugate reflection could be reduced to the 
extent that the signal is displaced from the probe beam path, diminishing 
the reflectivity measurement. 
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PROBE INTENSITY (W/cm2) 


Reflectivity variation with relative backward pump 
intensity. The dashed lines are fits to the three most 
intense probe beam data points for each backward pump 
intensity. 






















I 


An experiment was performed to maximize the transfer of forward pump 
energy into the phase conjugate return beam. In this experiment, the 
lenses used to focus the pump and probe beams in the four wave mixing cell 
were intercharged such that the probe beam completely enveloped the pumped 
volume in the cell. With this reconfiguration, the pump spot cross section 
was reduced to a 0.43 mm x 0.74 ran ellipse, while the probe cross section 
was increased to a 1.48 mm x 2.3 mm elipse. 

Figure 49 presents the variation of the four wave mixing efficiency 
with probe beam where the efficiency is defined as the fraction of the 
forward pump beam energy transferred into the signal beam. A maximum 
efficiency of 16% was demonstrated. 

The data is seen to follow the predicted efficiency, shown as a dashed 
line in the figure. The efficiency is given by 

I, 

V - 7-h- * (B-27) 

2 *3 

where the I 2 and I 3 are the backward pump and probe beam intensities, 
respectively. Note that the experiment was performed in the weak probe 
regime with no pump depletion. As with the reflectivity measurements, the 
experimental measurements outperformed the predictions in this regime. 

4.4.4 Conclusion 

Four wave stimulated Brillouin scattering phase conjugation has been 
demonstrated to given reflectivities as high as 100 , 000 % and efficiencies 
as high as 16%. The dependence of the reflectivity on pump intensities, 
relative backward pump intensity and probe intensity was measured and 
compared to simple theoretical predictions. 

The theory gave good estimates of reflectivity in the regime where I 3 
> I 2 . However, in two instances, the theory underestimated the 
reflectivity. First, the theory did not predict the large reflectivities 
we measured 1 3 / 1 2 « 1. In addition, the theory underestimated the 
increase in the reflectivity with pump intensity, as shown in Section 
4.4.2. In these experiments, the reflectivity was observed to grow with 
increasing 12 / 11 . 
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Figure 49. Average efficiency versus probe intensity. Efficiency 

was measured using a probe beam which enveloped the pumped 
volume in the four wave mixing cell. The theoretically 
predicted efficiency is given by the dashed line. 
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contrary to the prediction. .This disagreement between the theory and the 
experiment could be eliminated by the inclusion of backward pump intensity 
noise is the modelling. The experiment could recover the expected R a 
Ij/1 2 for I 2 well in excess of the background noise. 

Our attainment of reflectivities limited to lOOOx was probably due to 
our limited forward pump power. Higher pump powers would allow the invest 
igation of reflectivity scaling at giant reflectivites in the regime where 
R ~ Ij/I 2 , with I 2 far greater than the background roise. Investigations 
by Russian and British researches [23,36] have obtained much larger 
reflectivities (one-million fold) using much lower I 2 /Il with much larger 
forward pump Intensities (** 1 GW/cm^). 
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V. CONCLUSION 


In this program, three major concepts of phase conjugate PEL 
resonators were defined: the phase conjugate mirror, phase conjugate 
outcoupling and phase conjugate beam conversion. Each concept was 
developed to solve one or more problems that are pertinent to the 
conventional PEL resonator designs. Together, the new novel resonator 
designs are capable to satisfied the requirements for a ground-based high 
power PEL system. 


For each concept, a new novel unconventional adaptive optics technique 
is identified. The new nonlinear optics technology includes the effects of 
multiple-pulse stimulated Brillouin scattering, short-pulse four wave mix¬ 
ing and four-wave stimulated Brillouin scattering. Each technique was 
found to well satisfy the requirements of the corresponding resonator con¬ 
cept. 

Very few theoretical or experimental studies of these three effects 
have been reported before. A series of analytical theories, simulation 
codes and analyses are thus development to enhance our understanding of the 
effects and determine their feasibility to meet the system requirements. 


For verification and scaling purpose, proof-of-principle experiments 
have been conducted for the short-pulse four wave mixing and the four-wave 
stimulated Brillouin scattering phenomena using the picosecond and nano¬ 
second PEL simulators, respectively. A reflectivity of 10 3 has been 
observed for the latter phenomenon that demonstrated the feasibility of up- 
linking a high power conjugate beam from a ground-based FEL system. 

The research in this program provides a first step (a promising first 
step) toward the realization of a simple, reliable and powerful PEL system. 
The physical bases of the technology involved are mostly understood from 
the theoretical and experimental efforts under this program. The system 
concepts were evaluated and proven to be viable from a number of system 
, analyses. 


The research has been carried out by a team of excellent elite scient¬ 
ists in TRW. Therefore, it is not surprised that a number of invention 
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disclosures, publications, and conference presentations are associated 
exclusively with this program. We are proud to present the following list 

Invention Disclosures 

1. "Phase Conjugate Free Electron Laser Resonator", Docket number 
11-0440 

2. “Free Electron Laser Resonator With Phase Conjugate Outcoupling", 
Docket number 11-0439 

3. "Discriminative Scattering Scheme For Four Wave Mixing", Docket 
number 11-0438 

4. "Efficient Phase Conjugate Conversion of High Power Laser Beam", 
Docket number 11-0437 

Publications 

1. "Modeling of Transient Four-Wave Mixing", Proceeding Of SPIE Con¬ 
ference On Optics And Optoelectronic System. 

Volume 642, 1986. 

2. "Theory Of Short-Pulse Four Wave Mixing", submitted to Optics 
Letters, May, 1986 

3. “Analytical Theory Of Four Wave Stimulated Brillouin Scattering 
In The Saturation Regime", in process of clearance for publication 
in Physical Reviews. 

Conference Presentations 

1. "Modeling Of Transient Four-Wave Mixing", presented at SPIE Con¬ 
ference on Optics and Optoelectronic Systems, Orlando, 

April 2, 1986 

2. "Phase Conjugate Free Electron Laser Resonators", presented 

at International Quantum Electronics Conference, San Francisco, 

June 12, 1986. 

3. "Reflectivity Scaling With Fluence In Picosecond Four-Wave Mix¬ 
ing", presented at International Quantum Electronics 
Conference, San Francisco, June 13, 1986. 
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APPENDIX A 


Theory of Short-pulse Four Wave Mixing 

Chun-Ching Shih 
TRW, 01/1070 
One Space Park 
Redondo Beach, CA 90278 

Abstract 

An analytical solution of the transient four wave mixing phenomenon in a 
Kerr medium that is long compared to the pulse length is presented. 















V 



Optical phase conjugation via backward four wave mixing in a nonlinear 

medium can be described by two coupled-mode equations for the probe and the 

conjugated field.* 8y neglecting temporal derivatives in the equations, a 

steady-state solution has been found which gives the well-known reflectivity 

expression |tan(<L)| , where < is the field coupling strength and L is the 

2 

interaction length. This expression is the result of solving the coupled 
equations with time-independent boundary conditions and cannot be applied to 

3 

the transient four wave mixing involving time-varying pulses. 

Several attempts to solve the time-dependent coupled equations in the 

4-7 

transient regime have been reported in recent years. These solutions are 

based on the assumption of using either CW pump beams (constant k) or slowly- 

varying pump pulses that are much longer than the Interaction length 

(time-dependent <). Consequently, the result cannot be applied to the 

picosecond four wave mixing experiments that have received extensive attention 
8 9 

recently. ’ In these experiments, the probe and the pump pulses are supplied 
from the same picosecond laser source. The lengths of the interacting pulses 
are thus comparable with each other and, in many cases of interest, are 
shorter than the length of the nonlinear medium cell. The analysis of such 
experiments requires a fully transient theory of the short-pulse four wave 
mixing. 

In this letter, we present the first complete solution, ever reported in 
the literature, of the time-dependent coupled equations, the coupling 
strength, which is proportional to the product of the pump fields, is a 
function of position as well as of time. We found that, for most experimental 
situations in the laboratory, the reflectivity can be obtained in closed form, 
which is convenient for parametric studies of the phenomenon. The result can 
also be generalized to account for the bandwidth effect in the limit when the 
coherence length is much shorter than the pulse length. 

For convenience and simplicity, we will assume the slowly varying 
envelope approximation (SVEA) for the fields throughout this letter. We also 
assume weak probe and conjugated waves to avoid the complexity of pump 
depletion, self-focusing, and self-phase-modulation. Finally, we consider 
isotropic and instantaneous response of the Kerr effect in the nor-absorbing 
medium. These are the typical assumptions involved in the coupled-mode theory 

3 

of the four wave mixing effect. 

Consider the system consisting of four waves propagating in either 
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positive or negative z direction (Figure 1): 


A C (Z + vt) 


A 0 (Z - vt) 


Ai (Z - vt) 


,XV 2(Z+ 




Initial Conditions 

Ap(-°°. -°°) - Aq{Z - vt) 

A c (°°. -°°) - 0 


Final States 

Ap(°°, °°) - A p (Z - vt) 
A c {-°°, °°) - A C (Z + vt) 


Figure 1. Schematic of short-pulse four wave mixing. The pulse 
Interaction is assumed to be well within the medium 
cell so that the boundaries can be placed at z * 

• with Initial and final states defined at t - -• and 
•, respectively. 


EjU.t) 

E 2 (z.t) 

Ep(z.t) 


1/2 Aj(z - vt) exp[i«|t - k j z] + c.c., 
1/2 A 2 (z + vt) exp[iu 2 t + k 2 z] + c.c., 
1/2 Ap(z,t) exp[ioi p t - kpz] + c.c.. 


(A-l) 


E c (z,t) = 1/2 A c (z,t) exp[iu c t + k c z] + c.c.. 


where A^(z-vt) and A 2 (z+vt) represent the moving and non-depleting pump pulse 

profile, A and A are the slowly varying amplitude of the probe and the 
P C 

conjugated waves. By substituting Equation (1) into the inhomogeneous wave 
equation with nonlinear polarization and performing the following coordinate 
transformation: 


x = z - vt 


y = z + vt 


(A-2) 


















we obtain the time-dependent coupled equations for A p and A c : 


A c = 1eA 1 (x)A 2 (y)A p %xp[l(k p - kj)x + i(k 2 - k c )y] 
jy A* = i8*A*(x)A*(y)A c exp[i(k x - k p )x + l(k £ - k 2 )y], (A-3) 


where 


B-f jf x(3) (SI) 

« ^ X ^ (Gaussian). 


(A-4) 


The pump field amplitudes and A 2 , which are included in the constant k, are 
now written out explicitly. The exponential factors In Equation (3) represent 
the phase mismatch for the nondegenerate case. These factors, containing x 
and y, can be absorbed into the profiles, A^x) and A 2 (y), and do not affect 
the following derivations. For simplicity, the factors are dropped from now 
on. 

The coupled equations can be further simplified by defining the 
normalized pulse profiles. 


A c (x,y) = A c (x,y)/A 2 (y) 
A p (x ,y) = A p (x,y)/Aj(x). 


(A-5) 


and two fluence variables. 


x 

P = / 
00 


I Aj(x')|^dx' 


q = / lA 2 (y* )| 

- ao 



(A-6) 


In Equation (5), we note that the amplitude profile is normalized with respect 






















to only the pump pulse moving In the same direction. In Equation (6), p and q 
are the accumulated intensity (or fluence) from the front of the pump pulse to 
the point x or y defined by Equation (2). Since the interaction occurs well 
within the actual medium range, we can extend the medium boundaries to 
Infinity In both directions without affecting the result. Under this boundary 
condition, the initial and final values of p and q are. 


: p = o 
q = o, 


(A-7) 


: P * / |A,(x)r dx = P, 


q = / |A 2 (y)r dy = P ? , 


(A—8) 


where Pj and P^ are the pulse energy density of the pump beam 1 and 2. 

The equations now become 

3 ~ ~ it 

— A c (p,q) *-1 8 A (p,q) 


(A-9) 


— A s(p,q) = i B A (p,q), 
3q P 


with the initial conditions 


Ap(p,o) * A Q (p) 


A c (o,q) = o. 


(A—10) 


Equation (9) can be solved by using successive Laplace transforms with respect 
to p and q, 


[»■< 

















From Equation (11), we have 


W- S3^ ^> 

p q 




WV s TV 2 

P q 


(A-12) 


and finally. 



(A-13) 


where I„ and I, are the modified Bessel functions of zeroth and first order, 
o l . 

Equation (13) shows how the final pulse profiles can be calculated from the 
initial probe pulse profile. 

Consider the case that the probe and the codirectional pump pulses have 
the same profile except for a proportional constant, i.e.. 



(A-14) 


where Pp is the initial energy density of the probe pulse. This condition is 
most commonly seen in the laboratory when a single laser source is used. The 
final profile of the conjugated pulse can be readily found from Equation (13), 


A c ( p i«q) = i i fl | 


| 8 | 



I,(2|e|/PJq). 


(A-15) 


A (p. ,q) can be interpreted as the energy transfer rate from the beam 2 to 



















the conjugated beam. This rate is not uniform but depends on the fluence, q, 
of the beam 2. The nonuniformity can be seen by comparing the rate at the 
front and the rear end of the pulse. 


A c (Pi, p 2 )(Rear) Ij(2 |0| /P^) 

^ * 
A ( .(P 1 ,o)(Front) |6|/P^ 


(A-16) 


This ratio depends on the coupling strength and the product of pump energy 
densities. 

From Equation (15) we can find the reflectlvity, defined as the ratio of 
the energy densities in the conjugated and the probe pulses: 


P 00 

R = jr x ff- f dy | A.( -®,y) | 2 
P p -«• 

3 -p- f 2(ic iIA c ( p i»q)l 2 

P 0 

■ ! o (n) + (A-17) 

where n = 2 je| (A-18) 


From the result in Equation (17), we found that the reflectivity depends 
directly on the product of pump pulse energy densities, but not on the 
detailed structure of any pulses, as long as the pulses in the probe and the 
pump beam 1 have the same profiles. Choose a nonlinear medium, e.g. CS^» The 
dependent variable, n, for the case of equal pumps can be calculated as 

n (CS 2 ) = 63 Pj 2 (Joules/cm 2 ). (A-19) 

The dependence of R on this variable is shown in Figures 2 and 3. It has a 
quadratic dependence at small n and exponential dependence at large n. 


















Furthermore, no oscillation condition can be found at any level of energy 
density In contrast to the solution of CW four wave mixing. In Figure (3), 
also show the good agreement between this analytical result (solid line) and 
the simulation (dots) using the two-dimensional four wave mixing code 
developed at TRW.^ 
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Figure 2. Reflectivity verse n in linear scale (Equation (17)). 
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Figure 3. Reflectivity verse n Is logarithm scale (Equation (17)). 

The dots are the simulation result from a two- 
dimensional four wave mixing code (Reference 10). 

So far, we have neglected bandwidth effects and assumed only 
Fourier-transform limited pulses. It is of interest to generalize the 
solution to the cases where the coherence length is shorter than the pulse 
length. For this purpose, we consider a square pulse profile and divide the 
pulse into slots. The slot is equal to the coherent length. The fields in 
different slots are completely uncorrelated. Since the result in Equation 





















(17) Is applied to arbitrarily short pulses, we can sum up the Interactions 
between slots of different pulses to obtain the total reflectivity, 

R(r) = j [Ig(rn) + if (rn) - 1], (A-20) 

where r is the ratio of the coherence length to the pulse length. The 
reduction in reflectivity due to this bandwidth effect is shown in Figure 4 
for r = 0.1, 0.2, 0.4, 0.6, 0.8, and 1.0. The decrease is substantial at high 
n. At small u, R(r) is proportional to the ratio r. 



0 1 2 3 4 5 6 7 

n 

Figure 4. Reflectivity verse n for different ratios of coherence 
length and pulse length (Equation (20)). 

















In conclusion, we have successfully solved the time-dependent coupled 
equations for short-pulse four wave mixing. The theory has no restriction on 
the pulse profile of all participating beams. By using double Laplace 
transforms and a series of variable transformations, we are able to obtain the 
conjugate reflectivity in closed form. The result is then generalized to 
Include the bandwidth effect. 

The author acknowledges many stimulating discussions with Drs. D. L. 
Bullock and J. Menders in TRW. 
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Analytical Theory of Four Wave Stimulated Brillouin 
Scattering in the Saturation Regime 

Chun-Ching Shih 
TRW, 01/1070 
One Space Park 
Redondo Beach, CA 90278 

Abstract 

An analytical solution of the enhanced stimulated Brillouin scattering 
arranged in four wave mixing configuration is presented. The theory includes 
the effect of pump depletion and is able to demonstrate the reflectivity, 
efficiency and field variations in the saturation regime. 







I. Introduction 

Four wave mixing (FWM) via third-order optical susceptibilities in a 

nonlinear medium has been recognized as an attractive optical phase 

conjugation technique for its simple configuration, low power requirement and 

high reflectivity.* Although early theory of FWM predicted an infinitely 

2 

large reflection coefficient under certain pump conditions, the highest 

reflectivity ever observed has not exceeded 100. This is due to its 

stringent requirements which are difficult to meet in most experiments. Later 

analyses showed that this reflective singularity is highly unstable with 

4 5 

respect to the change of pump Intensity. * A small amount of pump depletion 
can drive the Interaction out of the resonance regime and thus reduce the 
reflectivity significantly. 

Recently, conjugate reflectivities in the range of 10 4 to 10 6 have been 
observed with the technique involving the combination of two major optical 
phase conjugation concepts: four wave mixing configuration and stimulated 
Brillouin scattering (SBS) process. 6,7 Together, we shall call it four wave 
stimulated Brillouin scattering (FWSBS). In this effect, two beams of the 
same frequency, w, propagate toward each other (Figure 1). A probe beam of 


*3 



slightly higher frequency is incident at a small angle into the interaction | 

region. The frequency difference should be close to the acoustic wave 
frequency, fi, of the backward SBS effect. This arrangement is identical to 
the ordinary nondegenerate FWM. However, the interaction mechanism is 

completely different. The SBS interaction between the backward pump (beam 2) | 

and the probe signal (beam 3) creates an acoustic grating. This grating is 
further enhanced by the stimulated scattering of the forward pump (beam 1) 
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into the return field (beam 4). The result is a substantial energy flow from 
beam 1 to beam 4 which is a phase conjugation of the probe field. 

To date, the theoretical studies of this effect have been limited to the 

8 9 

constant pump approximation. * The results showed that the FWSBS 
reflectivity has no singularity in the normal operational regime. Therefore, 
the energy transfer from the pump beams can be as complete as possible without 
affecting the scattering process significantly. It was also found that the 
relectivity depends only on the ratio, but not the product as in FWM, of the 
pump intensities. The reflectivity can be brought to an extremely high level, 
say 10®, by lowering the intensity in one of the pump beams. In order to 
demonstrate that the FWSBS can be a high reflectivity as well as a high 
efficiency process, the effect has to be studied in the saturation regime. 

In this paper, we report an analytical solution of the four coupled 
equations describing the wave interaction in FWSBS. 

In Section II, we derive the coupled equations and discuss the physical 
significance of each term. In Section III, we solve the equations with 
constant-pump approximation and check the conditions for high reflectivity. 

In Section IV, we give the analytical solutions of the four equations 
including the pump depletion effect and show the saturation phenomenon. 
Finally, a conclusion is drawn in Section V. 

11. Coupled Equations 

Consider the system consisting of four optical waves propagating in 
either positive or negative z direction and an acoustic wave with frequency, 


= E 1 exp[ia>jt + ikjz] + c.c. 

$2 - E 2 ex P[ lu 2 t: “ lk 2 Z ^ + C,C * 

= Ej exp[ioj^t + ikjz] + c.c. 

£4 = E 4 exp[iu> 4 t - ik 4 z] + c.c., 


(B-l) 


^ + kp + k3 + k 4 

d, = Q expfint + i - 4 - z] + c.c., 


(B-2) 


where Ej and Q are the slowly varying amplitude and the wave number of the 
acoustic wave is taken to be the average wave number of the four optical 
waves. For simplicity, the frequency difference between the pump beams and 














other beams is assumed to be in resonance with the phonon field. 


fl = U)-(U^ =0)3-0), 


(B-3) 


where 


(o = Wj = a>2. (B-4) 

Using Equations (1) and (2), the basic equations for the slowly varying 
amplitudes can be written as ^ 

(7 3 ? ' af) E l = ~ 1 9 1 Q * E 4 exp [-1 z] 

(7 af + 4 )E 2 = E 3 ex P I* 1 T- (B-5) 

^7 77 + 77^4 = E i ex P[i TT 

(al + r)Q = ig 2 {E^E 4 exp[-i $r- z] + E 2 E* expfi z]}, (B-6) 


where v is the light velocity in medium, and g 2 are the coupling 
coefficients, r is the phonon relaxation rate, and the phase-mismatch, 4k, is 
(Figure 2) 




















Figure 2. k-vector mismatch for the four waves 


In Equations (5) and (6), we have assumed that the optical waves are coupled 
to each other through the same phonon field with slightly different 
phase-mismatch conditions. For steady-state analyses, the time derivatives 
can be neglected and the equations become 


ir ■! >4 e i + ! hhh «pH* k *] 

IF = t ¥3 + f E l E 2 E 4 ex Pl l4kz l 

3E 

IT s f ! 1 E 4 + f E 1 E 2 E 3 ex P[ iAkz l 

where 

>j ■ |E / 

29^2 
9 = f * 


(B-7) 


(B-8) 


The two terms on the right hand side of Equation (7) have similar 

















V % 



expressions with the SBS and the nondegenerate FWM. The Interaction process 

of FWSBS can be explained from Equation (7). Assume the power of beam 1 is 

much higher than the other two input beams, but is still below the threshold 

of the ordinary SBS effect. Although beam 1 does not have sufficient power to 

induce a Stokes wave from the noise, it can be reflected from the acoustic 

grating, created already by the SBS interaction between beam 2 and beam 3. 

This part of the field is represented by the second term in Equation (7). The 

reflected wave accumulates along its path in the opposite direction of beam 3. 

At low interaction strengths, this parametric superposition process will go 

(3) 

through the whole medium, which is similar to the FWM process in a x medium 

except that the coupling coefficient is real, not imaginary. 

However, if the reflected wave becomes strong enough, the acoustic 

grating will be enhanced by the SBS interaction between beam 1 and beam 4. 

The reflected field is thus further amplified like in an ordinary SBS process 

that is represented by the first term in Equation (7). Due to the transition 

between different interaction processes, FWSBS has the combined advantages of 

the FWM and SBS effects: good conjugation fidelity because the reflected wave 

* 

is generated from the field mixing (E^ ~ E^E-j), not from noises; and high 
scattering efficiency because the amplification is a stimulated process that 
avoids the phase-mismatch problem. However, it is noted that FWSBS is a 
threshold phenomenon because of the stimulated amplification. For the 
conjugation process to be dominant, the second term in Equation (7) should be 
larger than the first term where E^ starts from the noise level. This 
condition defines the threshold level for 1^, 

l iU 

I 3 > ^"- (B-9) 

Consider the noise level is roughly thirteen orders of magnitude lower than 
the pump intensity. The threshold condition for the three input beams becomes 

V 3 > 10 ‘ 13 >!• (B-10) 

III Constant-Pump Approximation 

If the probe intensity, 1^, is much lower than the pump intensities, 1^ 
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and I 2 , the first two equations in (7) can be neglected. The remaining two 
coupled equations for Ej and E 4 can be solved exactly assuming constant Ij and 
Ig.®*® Consider the nonlinear medium extends from 0 to L in the z direction. 
Using the boundary condition, E 4 (0) * 0, the reflectivity was found 


I 4 U) 

TvtfT 


|o - 8 coth 6 | 


(B-ll) 


where 


J j = ® LI J 


2 . i i a/2 


0 - (a + JiJjO 


(B-12) 


J. - J 2 

o * ■— - » —- + iAkL. 


The dimensionless variable, J^, is a measure of the intensity, Ij, scaled by a 
factor gL. For a typical medium, acetone, at 1pm wavelength, we calculate Ak 
= 2 cm~*. In the limit of short interaction length, AkL << 1, the 
reflectivity can be further simplified. 


R = - 


I t I 2 [exp(J 1 + J 2 )-l] 2 
[I x + I 2 exp(Jj + J 2 )j‘ 


(B-13) 


It can be seen that the reflectivity is finite for all possible Ij and L 
and does not rely on unstable singularities. Now, from Equation (13), we can 
determine the condition of pump intensities to obtain high reflectivities. Tc 
avoid the ordinary S 8 S effect, it is required that Jj 2 < 30 * a threshold 
constant.** For example, we can choose = 20. At this power level, the 
exponential factor is exp(J 1 ) ~ 5 x 10 8 which is still much larger than the 
ratio (I./I 0 ) in most experimental conditions. Therefore, it is a very good 
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approximation to let the exponential factor approach infinity and we obtain 
the reflectivity 



Equation (14) shows that the constant-pump reflectivity is simply equal 
to the ratio of the two pump beam intensities. High reflectivity can thus be 
obtained by just lowering the backward pump intensity. This dependence is 
interesting comparing to the behavior in FWM, where the reflectivity increases 
with both pump intensities. We should also point out that the intensities 
have the following relation 


! 4 ! 3 

T[ s T 2 ' (B—15) 


which Is similar to the equality of reflectivities for the SBS effect of two 
pump beams sharing the same interaction region. From Equation (15), we can 
see that if Ig approaches l 2 , a substantial energy transfer from pump beam 1 
should result. In this pump depletion regime, the reflectivity will saturate 
with increasing probe intensity and the result in this section cannot apply. 

IV Saturation Phenomenon 

The saturation effect of FWSBS should be studied from the complete set of 
coupled equations in (7). To simplify the problem, we neglect the exponential 
factor of phase-mismatch. This factor is important only when the four wave 
mixing dominates the process for a long distance. In the saturation regime, 
the FWM-dominant region is confined to a small area near the rear end of the 
medium. Within this region, the exponential factor is nearly a constant and 
can be dropped from the equations. With this approximation, the coupled 
equations become 


3E 1 _ g g * 

3z ~ 7 4 E 1 7 E 3 E 4 E 2 


•N 
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Equation (16) has the following constants, 

P 8 'l ' >4 

q * h - r 3 
r « El E 2 - E 3 E 4 

S * E 1 E 3 ' E 2 E 4* 

These four constants are not independent. 

Ir | 2 - fs | 2 = pq. 


(B-17) 

They are related to each other by 

(B-18) 


The four coupled equations are thus reduced to three independent constants 
plus an equation for any of the four intensities, for example. 


3I 4 

IT = gI l I 4 + gE 1 E 2 E 3 E 4* 

Define the average intensities Ij 4 and 1 23 * 


(B-19) 
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J 14 * 7 ^1 + V 
! 23 “ 7 ^2 + * 3 )* 


(B-20) 


From Equations (17), (19) and (20), we found the equtions for and I 23 . 


TT ■ h « r2 * s2 * P 2 )('l4 - £ ♦ I 14 Clf 4 - #) 1/2 } 

P 

*-£■ - h k -- 2 ** 2 ♦ Aols - £ - iodfj - #) 1 / 2 f- 


(B-21) 


These two nonlinear equations can be further simplified by using the following 
transformation. 



(B-22) 


The equations for x and y are 

|f - Tjij- ( [ ( r * s) 2 ♦ p 2 ] x 2 - [(r - s) 2 * p 2 ]) 

([(r - s) 2 * q 2 ] , 2 - [(r * s) 2 ♦ , 2 ]). 


(B-23) 


These two equations can be integrated readily. The results are 










e' 


x + A * x~ + A 


y - B = y o ~ B yz 

y + B y + B * 
J J o 


(B-24) 


where 

a 2 - Cr - s) 2 + p 2 
(r ; s) 2 * P 2 


I 2 ■ ^ [p 2 ♦ q 2 * 2(p 2 ♦ s 2 )]. 


(B-25) 


and x Q and y Q are the values of x and y evaluated at z * 0. The solutions in 
Equation (24), x(z) and y(z), can be substituted into Equation (22) to obtain 
the average intensities, I^(z) and 123 ( z )* ^ he intensity of each beam can 
thus be calculated 





















Equations (22), (24) and (26) represent the complete analytical solution 
for the FWSBS effect. The result is expressed in terms of the boundary 
condition at one end of the medium. In reality, however, the boundary values 
are Ij(L), I 2 (0), I 3 (L) and I 4 (0), which are defined at either z = 0 or z = L. 
It is very difficult to express the constants, p, q, r, s, in terms of these 
known quantities. 

For the purpose of demonstration, we write Equations (22), (24) and (26) 
in a code which calculates 1.(0) from the input values I -(L). By fixing ML) 

J J • 

and I 3 (L), we can adjust the inputs I 2 (L) and I 4 (L) such that the output 
values of beam 2 and beam 4 equal the boundary condition I 2 (0) and I 4 (0). 

This process is quite stable and converges very quickly. 

To demonstrate the saturation effect, we choose the forward pump 

2 2 
intensity to be 500 MW/cm and the backward pump intensity 0.5 MW/cm . In a 2 

cm long acetone medium (g = 0.02 cm/MW), these correspond to Jj = 20 and J 2 = 

0.02, respectively. The probe intensity is varied over three orders of 

magnitude from 2.5 kw/cm to 2.5 MW/cm (Figure 3). As anticipated, the 
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Figure 3. Saturation effect of FWSBS. The curves are found 
to be well described by R = I./(I ? + I,) and 
n = Io/(I ? + U). 
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reflectivity in the non-depletion regime approaches 10 which is the ratio of 
the pump Intensities, yy The saturation Is observed as the reflectivity 
decreases from 90 percent to 10 percent in the range of I 3 = 0.1 I 2 to I 3 = 10 
I 2# The dependence seems to be well described by the following formulas 


R 



1 

+ I 
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n 


I 


3 




(B-27) 


where the efficiency, n, is defined as the fraction of the forward pump 
intensity transferred to the reflected wave. 

Equation (27) can also be derived without actually solving the coupled 
equations. Since s is a constant and E^(0) = 0, the fields at both boundaries 
have the relation 

E l (0)E 3 (0) - E 1 (L)E 3 (L) - e 2 (l)e 4 <l). (b _ 28) 

In most cases of interest, E 3 (0) is several orders of magnitude lower than 
other quantities and the lefthand side of Equation (28) can be neglected. The 
reflectivity is then found to be 


yu yu 
R 5 Tjfuy = T~rry* 


(B-29) 


Since 13 ( 0 ) is very small, I 2 (L) is approximately equal to [I 2 (0) + I 3 O-)] 
which gives the formula in Equation (27). Again, Equation (29) is an example 
of equal reflectivities discussed in Reference 12, where E 3 (0) is at the noise 
level. 

Equation (27) shows that high reflectivity and high efficiency can be 

obtained by properly choosing the input intensities. For example, if and 

Io are lower by two orders of magnitude from the case in Figure 3, the 
J 5 

reflectivity will be approaching 10 and the efficiency is still the same. At 
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these intensity levels, the phase conjugation threshold condition in Equation 
(10) Is still satisfied. 

The intensity variation in the medium is demonstrated for the case of I 2 
= Ij where the reflectivity is half of Its non-depletion value and the 
efficiency is near 50 percent (Figure 4). It is observed that Ij(z) and ^(z) 



Figure 4. Intensity variation in the interaction region for the case of ^ * 
Ij in figure 3. The four wave mix process dominates only in the 
first one-fifth of the region. After that, the reflected wave 
grpwm exponentially. 














have similar dependences except for a constant factor, this can be understood 
because x(z) and y(z) in Equation (24) have an identical structure. For most 
of the interaction region, the reflected wave intensity, I^(z), is subject to 
an exponential growth. This is the SBS-dominant region. Near the rear end of 
the interaction region (about one-fifth of the total length), the FWM process 
dominates the field growth. This further justifies the drop of the 
phase-mismatch exponential factor at the beginning of this section because the 
actual interaction length for FWM is much shorter than the medium length. 

So far, we have assumed that the optical waves are coupled to each other 
through the phonon field only and neglected the involvement of other nonlinear 
effects, such as Since the second term in Equation (7) is similar to 

the structure of the coupled equations for FWM in a Kerr medium, we expect 
that the coupling through x^ should take place. For completeness, the 
coefficient of the second term in Equation (7) should be replaced by 

g - g ± 1 8 n 2 k, (B _ 30) 

where n 2 is defined in the effective index of refraction 

" ■ "o + "2 <’>• (B-31) 

and <I> is the time average of the total intensity. In Equation (30), we add 
an imaginary part to the coupling coefficient. The positive sign is for beams 
1 and 3, and negative sign for beams 2 and 4. Consider the medium CS 9 , where 

O p ^ 

n 2 is 2 x 10 (cm /MW) and g is 0.2 (cm/MW). We found that the imaginary 
part is one order of magnitude smaller than the real part in Equation (30). 
Therefore, it seems appropriate to neglect the contribution of FWM through the 
third-order nonlinear susceptibility in the FWSBS effect. 

V. Conclusion 

We have derived an analytical theory for the four wave stimulated 
Brillouin scattering effect in the saturation regime. The analytical solution 
shows a moderate saturation effect due to the pump depletion. Because of its 
high conjugation reflectivity and good energy transfer efficiency, the FWSBS 
is of interest to many applications which are the subject of further 
investigation. 


137 

















The author appreciates many delightful discussions with Dr. J. Menders. 

This research is supported by the Office of Naval Research, Pasadena. 

References 

1. D. M. Pepper and A. Yariv, "Optical phase conjugation using three-wave 
and four-wave mixing via elastic photon scattering in transparent media" 
in Optical Phase Conjugation , edi. by R. A. Fisher, Chapter 2, Academic 
press, N.Y. 1983. 

2. A. Yariv and D. M. Pepper, "Amplified reflection, phase conjugation and 
oscillation in degenerate four-wave mixing". Opt. Lett. 1_, 16 (1977). 

3. D. M. Bloom, P. T. Liao and N. P. Economou, "Observation of amplified 
reflection by degenerate four-wave mixing in atomic sodium vapor". Opt. 
Lett., 2, 58 (1978). 

4. W. P. Brown, "Absorption and depletion effects on degenerate four-wave 
mixing in homogeneously broadened absorbers", J. Opt. Soc. Am., 73 , 
629-634 (1983). 

5. M. T, Gruneisen, A. L. Gaeta, and R. W. Boyd, "Exact theory of pump 
wave propagation and its effect on degenerate four-wave mixing in 
saturable-absorbing media", J. Opt. Soc. Am., Bj?, 1117-1121 (1985). 

6. N. F. Andreev, V. I. Bespalov, A. M. Kiselev, A. Z. Matveev, G. A. 
Pasmanik and A. A. Shilov, "Wave-front inversion of weak optical signals 
with a large reflection coefficient", JETP Lett., 32, 625-629 (1980). 

7. A. M. Scott and M. Hazel 1, "Efficient phase conjugation by Brillouin 
enhanced four-wave mixing", OSA annual meeting, Washington, D.C., 

October 17, 1985. 

8. N. F. Andreev, V. I. Bespalov, A. M. Kiselev, G. A. Pasmanik and 

A. A. Shilov, "Raman interaction in the field of opposing light waves", 
Sov. Phys. JETP, 55, 612-617 (1982). 

9. A. M. Scott, "Efficient Phase conjugation by Brillouin enhanced four 
wave mixing". Opt. Commun., 4£, 127-132 (1983). 

10. G. C. Valley, "A review of stimulated Brillouin scattering excited 
with a broad-band pump laser", IEEE J. Quant. Elect., QE-22 , 704-712 
(1986). 

11. I. D. Carr and D. C. Hanna, "Performance of a Nd:YAG oscillator/ 
amplifier with phase conjugation via stimulated Brillouin scattering", 
Appl. Phys., B36, 83-92 (1985). 












12. R. L. Abrams, C. R. Giuliano and 0. F. Lam, "On the equality of 
stimulated Brillouin scattering reflectivity to conjugate reflectivity 
of a weak probe beam". Opt. Lett., 6 , 131-132 (1981). 

13. A Hardy and Y. Silberberg, "Saturation effects in phase-conjugate 
lasers", J. Opt. Soc. Am., 73, 594-599 (1983). 























